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According to cosmological inflation, the inhomogeneities in our universe are of quantum mechani- 
cal origin. This scenario is phenomenologically very appealing as it solves the puzzles of the standard 
hot big bang model and naturally explains why the spectrum of cosmological perturbations is almost 
scale invariant. It is also an ideal playground to discuss deep questions among which is the quantum 
measurement problem in a cosmological context. Although the large squeezing of the quantum state 
of the perturbations and the phenomenon of decoherence explain many aspects of the quantum to 
classical transition, it remains to understand how a specific outcome can be produced in the early 
universe, in the absence of any observer. The Continuous Spontaneous Localization (CSL) approach 
to quantum mechanics attempts to solve the quantum measurement question in a general context. 
In this framework, the wavefunction collapse is caused by adding new non linear and stochastic 
terms to the Schrodinger equation. In this paper, we apply this theory to inflation, which amounts 
to solving the CSL parametric oscillator case. We choose the wavefunction collapse to occur on 
an eigenstate of the Mukhanov-Sasaki variable and discuss the corresponding modified Schrodinger 
equation. Then, we compute the power spectrum of the perturbations and show that it acquires 
a universal shape with two branches, one which remains scale invariant and one with — i, a, 
spectral index in obvious contradiction with the Cosmic Microwave Background (CMB) anisotropy 
observations. The requirement that the non-scale invariant part be outside the observational win- 
dow puts stringent constraints on the parameter controlling the deviations from ordinary quantum 
mechanics. Due to the absence of a CSL amplification mechanism in field theory, this has also for 
consequence that the collapse mechanism of the inflationary fluctuations is not efficient. Then, we 
determine the collapse time. On small scales the collapse is almost instantaneous, and we recover 
exactly the behavior of the CSL harmonic oscillator (a case for which we present new results), 
whereas, on large scales, we find that the collapse is delayed and can take several e-folds to happen. 
We conclude that recovering the observational successes of infiation and, at the same time, reaching 
a satisfactory resolution of the inflationary "macro-objectification" issue seems problematic in the 
framework considered here. This work also provides a complete solution to the CSL parametric 
oscillator system, a topic we suggest could play a very important role to further constrain the CSL 
parameters. Our results illustrate the remarkable power of infiation and cosmology to constrain new 
physics. 

PACS numbers: 98.80.Cq, 98.80.Qc, 03.65.Ta, 03.65.Yz, 05.40.-a 



I. INTRODUCTION 

Inflation is currently the leading paradigm for explain- 
ing the physical conditions that prevailed in the very 
early universe [1-5]. It solves the puzzles of the stan- 
dard hot big bang phase and it explains the origin of 
the inhomogeneities in our universe [6-11] (for reviews, 
see Refs. [12-18]). According to the inflationary scenario, 
these inhomogeneities result from the amplification of the 
unavoidable vacuum quantum fluctuations of the gravi- 
tational and inflaton fields during a phase of accelerated 
expansion. In particular, inflation predicts an almost 
scale invariant power spectrum for the cosmological fluc- 
tuations [19], a prediction which fits very well the high 
accuracy astrophysical data now at our disposal [20-2G]. 

Often less emphasized is the fact that inflation is also 
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particularly remarkable from the theoretical point of 
view. Indeed, the inflationary mechanism for the produc- 
tion of cosmological perturbations makes use of general 
relativity and quantum mechanics, two theories that are 
notoriously difficult to combine. Moreover, this mecha- 
nism leads to theoretical predictions that are possible to 
study observationally with great accuracy. In fact, infla- 
tion is probably the only case in physics where an effect 
based on general relativity and quantum mechanics leads 
to predictions that, given our present day technological 
capabilities, can be tested experimentally. 

The situation described above can be used to inves- 
tigate deep questions. Among these deep questions is 
how the quantum measurement problem looks in a cos- 
mological context. According to inflation, the Cosmic 
Microwave Background (CMB) radiation anisotropy [27] 
is an observable and is therefore described by a quan- 
tum operator. As a consequence, when one looks at a 
CMB map, one observes the result of a measurement of 
that observable. According to the postulates of quan- 
tum mechanics in the Copenhagen interpretation, this 
means that the wavefunction of the inflationary pertur- 
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bations has collapsed to an eigenvector of this operator 
and that the CMB map corresponds to one of its eigen- 
value. The problem with this approach is that the col- 
lapse is supposed to occur only when an observer per- 
forms a measurement on the system. Clearly, there was 
no observer before or when the CMB was emitted. This 
seems to contradict the phenomenological fact that large- 
scale structure formation started early in the history of 
the universe since these structures are seeded by the same 
early physics which led to CMB fluctuations. As a mat- 
ter of fact, CMB fluctuations can also be understood as 
the earliest hint that primordial inhomogeneities had al- 
ready started to grow at that time. Furthermore, in some 
sense, the observers are actually the end product of the 
structure formation process! Of course, this measure- 
ment problem is already present in conventional labora- 
tory situations but it seems to be exacerbated (to use the 
words of Ref. [28]) in a cosmological context. 

Important steps towards a better understanding of 
these issues have already been accomplished. In par- 
ticular, it was shown that the inflationary accelerated 
expansion transforms a coherent vacuum state into a 
strongly squeezed state [29], the corresponding squeez- 
ing being much more important than whatever can be 
realized in the laboratory [30]. In this limit, the pre- 
dictions of the quantum formalism are indistinguishable 
from that of a theory where the fluctuations are just 
assumed to be realizations of a classical stochastic pro- 
cess [31-33]. The classical limit is a subtle concept in 
quantum mechanics but, in this sense (and in this sense 
only!), the system can be characterized as being clas- 
sical [34]. Moreover, the large-scale cosmological per- 
turbations are not isolated and, as a consequence, the 
phenomenon of decoherence [35-37] is relevant for them. 
This has for consequence that their density matrix be- 
comes diagonal before recombination, a criterion which 
is also considered as necessary in order to understand 
the quantum-to-classical transition [31, 32, 38-43]. How- 
ever, it is known that decoherence per se does not solve 
the measurement problem [44, 45]. Indeed, it remains 
to understand how a single outcome can be produced. 
This point is particularly important given that we only 
have one CMB map, that is to say only one measure- 
ment of the corresponding observable. In other words, 
even if the cosmological fluctuations can be viewed as a 
classical stochastic problem, this does not explain how a 
given realization of this process becomes an actual per- 
ception. This "macro-objectivation" problem is already 
present in a conventional situation but, as already men- 
tioned before, it becomes particularly embarrassing in 
the context of inflation where the collapse of the wave- 
function cannot be due to the presence of a conscious 
observer. Facing this situation, the common attitude is 
to postulate that decoherence should be combined with 
a new interpretational scheme, different from the Copen- 
hagen interpretation. Typically, in cosmology, the many 
world approach is often implicitly assumed [34, 4G-5(J]. 
Another frequently mentioned possibility, which seems to 



be particularly well suited to the cosmological context, is 
to consider that the wavefunction only represents the in- 
formation that we have on the system [51]. In this case, 
the issue of the wavefunction collapse becomes irrelevant 
since it just corresponds to a situation where the observer 
updates its knowledge (in the Bayesian sense) about the 
physical properties of the system. Other attempts, such 
as the non-local hidden variable theories, have also been 
tried [52-57]. In all these cases, the cosmological situa- 
tion does not differ much from a conventional laboratory 
situation and, moreover, does not lead to new, falsifiable, 
predictions^. Then, it becomes a question of taste which 
approach best fits one's own prejudices. 

However, there exists an exception to the conclusion of 
the previous discussion, namely the case of the collapse 
models [Gl-GG] (for reviews, see Refs. [G7, G8]). In this ap- 
proach, the Schrodinger equation is modified by adding 
non linear and stochastic terms which renders dynamical 
the collapse of the wavefunction. The model has nice fea- 
tures: firstly, the approach seems to follow a conservative 
strategy since, in physics, it is standard to first consider a 
linear theory and then, in order to have a more accurate 
description, to consider non linear corrections; in some 
sense, the collapse theories follow this line of argument. 
Secondly, there is now a single law of evolution for the 
state vector and, thirdly, the Born laws can be derived 
instead of postulated. There are also disadvantages such 
as the property that energy is not conserved or the fact 
that the relativistic formulation of the theory appears to 
be technically and conceptually difficult to develop (how- 
ever, see Ref. [69]). But, clearly, the main advantage in 
comparison to the possibilities discussed above is that 
this approach is falsifiable since it leads to predictions 
different from that of conventional quantum mechanics. 
This fact has been widely used in order to constrain col- 
lapse theories in the laboratory [G8, 70-73] but, clearly, 
it is also important to see whether this could be done in 
a cosmological context [74-7G]. It is therefore interest- 
ing to investigate what the collapse theories have to say 
about the infiationary mechanism. Notice that, regard- 
less of one's opinion about collapse theories, the subject is 
worth studying: a supporter would argue that the cosmo- 
logical measurement problem can possibly find a natural 
solution within this theory and an opponent would hope 
that the constraints obtained in a cosmological context 
can rule out the theory. In fact, this last question turns 
out to be very important. Indeed, as already mentioned, 
the constraints that exist on collapse theories are usually 
obtained from physical phenomena that can be observed 
in the laboratory. Therefore, by studying collapse the- 
ories in the context of cosmology and inffation, one can 
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hope to derive very relevant new constraints since one 
now deals with characteristic scales (energy, length etc 
. . . ) which typically differ by many orders of magnitude 
from those used to in a down to earth context. This il- 
lustrates again the conceptual relevance of inflation when 
it comes to very fundamental questions and its power to 
constrain alternatives to gravity but also to quantum me- 
chanics. In some sense, inflation represents an ideal play- 
ground to test new theories. Notice in passing that the 
very same strategy was used in the case of the so-called 
trans-Planckian problem of inflation [77-79] where it was 
shown that the inflationary observables could possibly 
contain an imprint (although probably small) of string 
theory. 

We are using (modified) Schrodinger-type of equation 
to describe the behavior of cosmological perturbations. 
This is justified because each Fourier mode of those ef- 
fectively evolves in an independant way and cosmologi- 
cal expansion permits to define a priviledged time. This 
allows for a sensible treatment of cosmological perturba- 
tions even though a fully relativistic CSL model, which 
could be naively expected to be required, is still lacking. 
At this moment, surprizingly, it is easier to treat infla- 
tionary perturbations than ordinary particle physics. 

It should also be emphasized that the idea of ap- 
plying collapse theories to inflationary perturbations 
of quantum-mechanical origin was first considered in 
Refs. [80-82]. In these articles, a phenomenological 
model for the collapse process was assumed and the corre- 
sponding physical properties were derived. In particular, 
the power spectrum of the perturbations was calculated 
and was shown to deviate from the standard predictions. 
Therefore, Refs. [80-82] have demonstrated that, in prin- 
ciple, it is possible to observationally test collapse theo- 
ries in a cosmological context. Our approach differs from 
that of Refs. [80-82] in the fact that we use the Continu- 
ous Spontaneous Localization (CSL) model to implement 
the collapse dynamics. This has the advantage that our 
calculations can be directly confronted and compared to 
other results obtained in other branches of physics. 

This paper is organized as follows. In the next sec- 
tion. Sec. II, we present a brief review of the theory of 
inflationary cosmological perturbations of quantum me- 
chanical origin. We especially focus on the calculation of 
the power spectrum since this quantity is the tool that 
allows us to relate the inflationary theory with the CMB 
observations. Then, in Sec. Ill, we discuss the cosmo- 
logical measurement problem and we explain how high 
accuracy CMB measurements can constrain inflation. In 
Sec. IV, we consider collapse theories, in particular, its 
CSL version, which is, as already mentioned, the case we 
use in this article. These sections aim at rendering the 
present work self-contained for readers with different ex- 
pertises. Then, we show how the harmonic oscillator can 
be treated in this context. This case is particularly rel- 
evant for cosmological fluctuations since it corresponds 
to the small-scale limit (in comparison to the Hubble 
radius) of the theory of cosmological perturbations. In 



Sec. V, we apply the CSL theory to inflation and to the 
calculation of the power spectrum. We use this result 
to constrain the parameter that controls the deviations 
from ordinary quantum mechanics. In Sec. VI, we study 
in more details the collapse phenomenon and explicitly 
compute the collapse time on small and large scales. In 
Sec. VII, we summarize our results and present our con- 
clusions. We end the paper with an Appendix A where 
it is shown that changing the "temporal gauge" in which 
the modified Schrodinger equation is written does not af- 
fect the shape of the power spectrum. This calculation 
reinforces the generic character of the results obtained in 
this work. 



II. INFLATIONARY COSMOLOGICAL 
PERTURBATIONS 

A. Basic Formalism 

By definition, inflation is a phase of accelerated ex- 
pansion that took place in the very early Universe, prior 
to the standard hot Big-Bang phase [1-5] (for reviews, 
see Refs. [13-15]). As is well known, postulating such 
a phase of evolution allows us to solve the standard 
problems of the hot Big-Bang model. Given that, at 
very high energies, field theory is the relevant frame- 
work to describe matter, a natural way to realize in- 
flation is to consider that a real scalar field (the "in- 
flaton" field) dominated the energy density budget of 
matter in the early Universe. Moreover, this assumption 
is compatible with the observed homogeneity, isotropy 
and flatness of the early Universe. Technically, the 
above-mentioned situation can be described by the met- 
ric tensor ds^ = — di^ + a'^{t)Sijdx^dx^ , where a{t) is 
the Friedman-Lemaitre-Robertson- Walker (FLRW) scale 
factor and t the cosmic time^. The Einstein equations 
imply that d/a = —{p+3p)/{6M^^), p andp being the en- 
ergy density and pressure of the matter sourcing the grav- 
itational field and M^^ the Planck mass (a dot denotes a 
derivative with respect to the cosmic time t) . For a scalar 
field, this reduces to &/a = V{if){l - (p'^ /V)/{3]VPJ, 
where V{(p) is the scalar field potential. This means that 
inflation (i.e. a > 0) can be obtained provided the infla- 
ton slowly rolls down its potential so that its potential 
energy dominates over its kinetic energy. This also shows 
that the inflaton potential must be sufficiently flat, a re- 
quirement which is not always easy to obtain in realis- 
tic situations and makes the inflationary model building 
problem a difficult issue [83]. The physical nature of 
the inflaton field has not been identified (there are many 
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candidates) and, as a consequence, the shape of V{(p) is 
not known. Of course, different V{(p) lead to different 
inflationary expansions but, since these different poten- 
tials must all be sufficiently flat, the corresponding scale 
factors are all approximately given by de Sitter solution. 
This solution is described by the scale factor a{t) ~ e^*, 
where H = a/a is the Hubble parameter, a slowly- varying 
quantity directly related to the energy scale of inflation. 
Observationally, this last quantity is not known but is 
constrained [22] to be between the Grand Unified The- 
ory (GUT) scale, that is to say ~ lO^^ GeV, and - 1 TeV. 
The previous considerations show that inflation can also 
be viewed as a phase of quasi-exponential expansion. 

A concrete illustration of the above discussion con- 
sists in considering power-law inflation [<S4]. Although 
it is based on a specific model with potential V{ip) — 
^./4g-Qv/J>^Pi (yfHYi Q. constant), it captures, in a simple 
way, all the essential properties of inflation and, more- 
over, is the only scenario which permits an exact inte- 
gration of the equations of motion (at the background 
level but also at the pcrturbative level, see below). The 
corresponding scale factor is given by 



1+/3 



(1) 



where £o is a length the value of which is fixed once the 
energy scale of inflation is known and rj in the conformal 
time defined by dt = ad?7, see Eq. (2). The quantity /3 
is a free parameter such that /3 < — 2 and is related to 
a through = (fi + 2)/(/3 -I- 1). The case /? = -2 

represents the de Sitter solution since it implies a = 0, 
i.e. a flat potential (and, of course, in cosmic time, the 
solution a cx 1/77 is given by an exponential). Therefore, 
different /? represents different inflationary solutions and 
/? must always be close to —2 in order for the potential 
to be sufficiently fiat. As announced, power-law infiation 
illustrates well the discussion of the previous paragraph. 

The above arguments can be considered as strong hints 
in favor of infiation. However, soon after its advent, it 
was realized that inflation, combined with quantum me- 
chanics, leads to an even more impressive result, namely 
it naturally explains the origin of the Cosmic Microwave 
Background (CMB) anisotropics and of the large-scale 
structures. According to the inflationary paradigm, these 
deviations from homogeneity and isotropy originate from 
the unavoidable zero-point quantum fluctuations of the 
coupled inflaton and gravitational flclds. Statistically, 
the fluctuations are characterized by their two-point cor- 
relation function or power spectrum. The observations 
[20-2G] indicate that the corresponding power spectrum 
is close to the Harrison-Zel'dovich, scale invariant, power 
spectrum with equal power on all scales. That this power 
spectrum represents a good flt to the astrophysical data 
was in fact realized before the advent of inflation but no 
convincing fundamental theory was known to explain this 
result. 

The main success of inflation is that it precisely pre- 
dicts an almost scale invariant power spectrum, the small 
deviations from scale invariance being connected with the 



micro-physics of inflation [6-1 1] . The fact that different 
types of inflationary scenarios lead to a power spectrum 
which is, at leading order, always close to scale invariance 
is connected with the fact that the inflationary scale fac- 
tor is always close to the de Sitter solution (see above) 
or, equivalently, with the fact that the inflaton potential 
is always almost flat. The deviations from scale invari- 
ance are related to the deviations from a flat potential 
and, therefore, depend on the detailed shape of the po- 
tential. As a consequence, measuring them allows us to 
say something about V{ip) and there is currently an im- 
portant effort in this direction using the high accuracy 
CMB data that have been released in the past years. 

Let us now see how the results reviewed before can 
be derived. Clearly, in order to model the cosmological 
fluctuations, one needs to go beyond homogeneity and 
isotropy. The most general metric describing small fluc- 
tuations of the scalar type on top of a FLRW Universe 
can be written as [12] 

ds^ = (ri) {- (1 - 2(/)) dr]^ + 2 {d^B) dx'dr] 

+ [(1 - 2tl;) 6,j + 2d,djE] dx'dx^ } . (2) 

A similar approach could be used to take into account 
tensor perturbations (i.e. gravity waves). Here, we do not 
include them since they are subdominant in the CMB, 
representing less than ^ 20 % at 2a confidence level [22] 
and, in addition, doing it would not bring any new as- 
pects to the question we want to investigate in this ar- 
ticle. In Eq. (2), the four functions (j), B, ip and E 
are of course functions of time and space since we con- 
sider an inhomogeneous and anisotropic situation. As 
is well known, the above approach is redundant because 
of gauge freedom [12, 85, 86]. A careful study of this 
question shows that the gravitational sector can in fact 
be described by a single, gauge-invariant, quantity, the 
Bardeen potential defined by [85] 



$3 {r,,x) = ^+-[a{B-E')]' , 



(3) 



where a prime denotes a derivative with respect to the 
conformal time 1]. In the same manner, the matter sector 
can be modeled by the gauge invariant fiuctuation of the 
scalar field 

(5v3(si) (^^ x) ^5ip + if' (B - E') . (4) 

The two quantities $3 and Sip'^^^^ are related by a per- 
turbed Einstein constraint. This implies that the scalar 
sector can in fact be described by a single quantity. For 
this reason, we now introduce the so-called Mukhanov- 
Sasaki variable [(i, 87] which is a combination of the 
Bardeen potential and of the gauge invariant field 



V {rj, x) = a 



6(p 



(gi) 



H 



(5) 



where H = a! j a. All the other relevant quantities can 
be expressed in terms of v{r\, x) which, therefore, fully 
characterizes the scalar sector. 
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The next step consists in deriving an equation of mo- 
tion for t)(?7, a;). This can be done directly from the per- 
turbed Einstein equations but, here, we first estabhsh the 
action for the quantity v{ri,x). Expanding the action of 
the system (i.e. Einstein-Hilbcrt action plus the action of 
a scalar field) up to second order in the perturbations, 
one obtains [f 2] 



(6) 

where ei = 1 — H' /H^ is the first slow-roll parame- 
ter [88, 89]. As the formula d/a = H^{1 — ei) shows, 
the condition ci < 1 is in fact sufhcicnt to have infla- 
tion. Moreover, we have slow-roll inflation [f9, 88-9f] 
if ei ^ 1. In this case, it is easy to show that ei ~ 
{M^j2V'^){dV/dip)'^, i.e. t\ is in fact a measure of how 
much the inflaton potential deviates from a flat potential. 
Equivalently, according to the previous considerations, 
this is also a measure of how much the inflationary expan- 
sion deviates from a pure de Sitter solution. In the case 
of power-law inflation, one has ei — (2 -I- /3)/(H- /3) and, 
of course, ei = when /3 = — 2 (de Sitter solution). The 
scale factor can also be rewritten as a(r\) ~ 
and this formula is in fact valid for any slow-roll model 
of inflation, i.e. for arbitrary shaped potentials, not nec- 
essarily of the exponential type. In this sense, power-law 
inflation with /3 < — 2 is a simple representative of all 
the slow-roll scenarios. Therefore, the fact that, in this 
paper, we focus on this particular model for technical 
reasons (again, because this model allows an easy inte- 
gration of the equations of motion at the background and 
perturbative level) does not restrict in any way the gen- 
erality of our considerations. 

Our next move consists in Fourier transforming the 
quantity v{r],x). This is of course justifled by the fact 
that we work with a linear theory and, hence, all the 
modes evolve independently. We have 



V {r],x) = 



1 



(27r 



,3/2 



(7) 



with v-k = "^fc because v{ri, x) is real. Then inserting 
this expansion into Eq. (6), one arrives at [12] 



(2)<55 = / d,j / d'k 



VkVk 



(8) 



where the integral over k is taken over half the Fourier 
space only. Next, we define pk, the variable canonically 
conjugate to 



Pk 



SC 



'^k > 



(9) 



where C is the Lagrangian density in Fourier space that 
can be derived from Eq. (8). This allows us to calculate 



the Hamiltonian which reads 



H 



d k i pkpl + Vkvl 



(10) 



This Hamiltonian represents a collection of paramet- 
ric oscillators (i.e. one oscillator per mode), the time- 
dependent frequency of which can be expressed as 



(11) 



We see that the frequency depends on the scale factors 
and its derivatives (up to the fourth). This means that 
different inflationary backgrounds (i.e. different inflaton 
potentials) lead to different w(ry, k) and, therefore, to dif- 
ferent behaviors for Vk{r]). From Eq. (10) or Eq. (8), it is 
easy to derive the equation of motion for the Mukhanov- 
Sasaki variable. One obtains 



+ uj^ (ri,k)vk = 0, 



(12) 



which confirms that each mode behaves as a parametric 
oscillator. Once a model of inflation has been chosen, the 
potential V{ip) is known and, hence, the corresponding 
scale factor can be calculated. This, in turn, allows us to 
determine uP'iji, k) and, then, one can solve the equation 
of motion (12). However, in order to flnd the solution for 
the Fourier component of the Mukhanov-Sasaki variable, 
one also needs to specify the initial conditions. Classi- 
cally, there does not seem to exist a natural criterion to 
choose them. However, when quantization has been per- 
formed, the requirement that it be initially in the vacuum 
state of the theory leads to well-defined initial conditions. 
We now turn to these questions. 



B. Quantization in the Schrodinger Picture 

In this section, we review how the cosmological pertur- 
bations are quantized. Very often in the literature, this 
is done in the Heisenberg picture. Here, we carry out the 
quantization in the Schrodinger picture [15] because this 
is more convenient for the problem we want to investigate 
in this article. In order to quantize the system, it is also 
more convenient to work with real variables. Therefore, 
we introduce the following definitions 



Vk 



1 

71 



Pk 



1 

71 



{pI + ^pD- (13) 



In the Schrodinger approach, the quantum state of the 
system is described by a wavefunctional, ^[v{ri,x)]. 
Since we work in Fourier space (and since the theory is 
still free in the sense that it does not contain terms with 
power higher than two in the Lagrangian) , the wavefunc- 
tional can also be factorized into mode components as 



^ Hrj, x)] ^ll^k {vl, vi) = n *^ (vl) (vi) 



(14) 



6 



Quantization is achieved by promoting Vk and to 
quantum operators, Vk and pk, and by requiring the 
canonical commutation relations 

[ik.Pq] - *<5 (fc - g) , [vlp'g] =i5{k-q). (15) 
These relations admit the following representation 

a* 



dv 



R,I 



(16) 



The wavcfunctional 5* [v{rj,x)] obeys the Schrodinger 
equation which, in this context, is a functional differ- 
ential equation. However, since each mode evolves in- 
dependently, this functional differential equation can be 
reduced to an infinite number of differential equations for 
each Concretely, we have 



R,I 



drj 



^ R, I 

where the Hamiltonian densities ' , are related to the 
Hamihonian by if = / d^fe {u^ + ii]^ . They can be 
expressed as 



(17) 



R,I 



1 



\^'{ri.k)(<^^')\ (18) 



where we have made use of the representations (16). 

We are now in a position where we can solve the 
Schrodinger equation. Let us consider the following 
Gaussian state 

(jl^ "fe'') = Nk{^)e-''''^^^^<'"y . (19) 

The functions Nkij]) and ^kiv) are time dependent and 
do not carry the subscripts "R" and/or "I" because they 
are the same for the wavefunctions of the real and imagi- 
nary parts of the Mukhanov-Sasaki variable (see below). 
Then, inserting 4'^. given by Eq. (19) into the Schrodinger 
equation (17) implies that Nk and flk obey the differen- 
tial equations 

TV' i 
ij^=^k, n',, = -2inl + -uj^{r^,k). (20) 



The solutions can be easily found and read 



Nk = 



. ^11. - - — ^ . 

2fk 



(21) 



where fk is a function obeying the equation f'k+uj'^fk = 
0, that is to say exactly Eq. (12). The first equation (21) 



guarantees that the wavefunction is properly normalized. 



^R,I^R,I*j^RJ ^ ^ (22) 

Let us now discuss the initial conditions. The funda- 
mental assumption of inflation is that the perturbations 
are initially in their ground state. At the beginning of in- 
flation, all the modes of astrophysical interest today have 
a physical wavelength smaller than the Hubble radius, 
i.e. k/{aH) — > oo. In this regime, one has ^^(77, k) — ^ fc^ 
and each mode now behaves as an harmonic oscillator (as 
opposed to a parametric oscillator in the generic case) 
with frequency lo — k. As a consequence, the differential 
equation for fkivi) can easily be solved and the solution 
reads fk = Ake^''^ + BfcC"*'^'', Ak and Bk being integra- 
tion constants. Upon using the second equation (21), one 
has 



2 AkC'-'''^ + SfcC-*'^') ■ 



(23) 



The wavefunction (19) represents the ground state wave- 
function of an harmonic oscillator if = fc/2. There- 
fore, one must choose the initial conditions such that 
Bk = 0. Moreover, it is easy to check that the Wronskian 
^ = f'kfk ~ fk fk is a conserved quantity, dW^/d?/ = 0, 
thanks to the equation of motion of fk ■ Straightforward 
calculation leads to = 2ik\Ak\^ ■ In the Heisenberg 
picture the canonical commutation relations require that 
W = i. Even if in the Schrodinger picture presently used, 
the specific value of W is irrelevant since it cancels out 
on all calculable physical quantities, this value is conven- 
tionally adopted, which amounts to setting Ak — l/\/2k. 
As announced, requiring the initial state to be the ground 
state has completely fixed the initial conditions. We 
see that Eq. (12) (or, equivalently, the equation for fk) 
should thus be solved with the boundary condition 



lim fk = 

k/(aH)^+oo 



1 



(24) 



This choice of initial conditions is referred to as the 
Bunch-Davies vacuum. 



C. The Power Spectrum 

Let us now turn to the calculation of the power spec- 
trum and first introduce the two-point correlation func- 
tion, defined by 



7 



(25) 



The next step consists in using the Fourier transform of the Mukhanov-Sasaki variable, see Eq. (7) and the exphcit 
form of the wavefunction of Eq. (19). One arrives at 

{^,\v{r^^x)v{n,x + r)\^)^-^ j dp dq e'^- e'-?' [] (^^^) / H ^^^'^ 5fen.[(.«)^+(.i 

I ""J J k ^ ^ k 

(26) 

If p 7^ ±<7, the result of the integration is zero since the integrand (up to the Gaussian weight) becomes linear in 
^R,i or v^'^. li p — q, then the only non-linear term in the integrand is given by {y^Y^ ~ (^p)^ Each term 
contributes the same amount, so the difference vanishes. The only possibility left is therefore p = —q, such that 
VpVq = (vpY + ("p)^ /2) the factor 1/2 coming from the definition of see Eqs. (13). This leads to 

{^\viv,xMri,x + rm = / dpe-^-n (^^) / Rd^.^d^e-^^^^'^^'^K^^^^+f-^^l «)^ (27) 

the factor of 2 originating from the fact that we have two contributions, one given by the term (wp)^ and the other 
by (fp)^. The Gaussian integrals can easily be carried out. They are all of the form "J da;e^"^ ", except of course 
the one over Vp which is of the form "J dx e~"^ " . As a consequence, one obtains 



\v(j], x)v{rj, x + r)\'i') 



(27r)3 



dpe 



-ip-r 



N 

n 



25Rel7fc\ 1 
2 

I 



N 

n 



N-1 

n 



(28) 



The infinite product "Hfe^ ^" means a product over all given scale, i.e. 
the wavevectors but p. One can always write this prod- 
uct as "•^23?e VLp/i: Hfc'" ; then the two last infinite prod- 
ucts in the above expression exactly cancel the first one. 
Therefore, we are left with 



n(fc) = ^IM 



(32) 



{^\v{ri,x)v{ri,x + r)\-^) = 



1 



(2^) 



dpe 



1 



43?e Q.p ■ 
(29) 

We now need to express SReOp in terms of the function 
/p. From the second Eq. (21), one easily shows that 



i W 



4|/pP' 



(30) 



and we obtain our final expression for the two-point cor- 
relation function 



(^E" \v{ri, x)v{ri, x + r)\ 4') = 



1 



(27r)3 



dpe-P--|/pp 



+ 00 



dp sinpr 
p pr 



(31) 



The same manipulations allow us to express the two- 
point correlation of two Fourier amplitudes. It can be 
written as 



\vkv*p\ *) = y n d<dz;^4'; Vk v*p *q 



(33) 



This integral is non-vanishing only if fe = p (otherwise 
one has to integrate an odd function) and receives two 

contributions, one from (v^)^ and the other from (wfc)^. 
Repeating calculations already performed before, one fi- 
nally arrives at 



27r2 



(34) 



We now need to explain how the cosmological pertur- 
where, in the last expression, we have used our choice bations of quantum-mechanical origin studied above are 
W = i. The power spectrum is just defined as the square related to observables in cosmology. This is the goal of 
of the Fourier amplitude per logarithmic interval at a the next section. 
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D. From Quantum Fluctuations to CMB 
Anisotropies 

The presence of quantum fluctuations in the infla- 
ton and gravitational fields has many observational im- 
plications. Here, we focus on one of them, namely 
the existence of CMB temperature anisotropies. The 
importance of this observable is that we now have at 
our disposal very high accuracy measurements of those 
anisotropies [20, 21]. Moreover, even more accurate data 
will be released soon [92]. The relation between the tem- 
perature fluctuations along a given direction e and the 
cosmological perturbations is expressed by the so-called 
Sachs- Wolfe effect [93, 94]. A simplified version of this 
result, valid on large angular scales only, can be written 
as [94] 



Vo) 



(35) 



where T represents the averaged background tempera- 
ture, i.e. T ~ 2.7 K, r^&s is the conformal time at emission 
(that is to say at the surface of last scattering) and ?7o 
is the present conformal time. The vector Xq landmarks 
the place of reception, in the present case Earth (or a 
satellite orbiting the Earth). The quantity C, denotes the 
curvature perturbation. It is related to the Bardeen po- 
tential defined in Eq. (3) through the following expression 
[12, 8G, 95] 



C 



2%- 
3~ 



1 



(36) 



where w = p/ p is the equation of state parameter, that 
is to say the energy density to pressure ratio of the dom- 
inant fluid. For instance, for the matter dominated era 
{w = 0), during which recombination takes place (at a 
redshift of z&s — 1100), on large scales, one simply has 
C ~ 5<I>3/3 since the Bardeen potential is constant. The 
importance of C, lies in the fact that it is a conserved 
quantity on large scales [86, 95]. Therefore, its spectrum, 
calculated at the end of inflation, can directly be prop- 
agated to the recombination time as it is not sensitive 
to the details of the cosmological evolution, in particular 
to those of the complicated reheating era [96-100]. The 
curvature perturbation can also be expressed in terms of 
the Mukhanov-Sasaki variable as 



c 



1 



(37) 



Finally, in the framework of the theory of inflationary 
cosmological perturbations of quantum-mechanical ori- 
gin, we have seen that v is in fact an operator. This im- 
plies that C and 5T /T are also quantum operators and, 
for this reason, from now on, wc will denote them with a 
hat. 

Since the operator ST jT lives on the celestial sphere, 
it can be expanded over the spherical harmonic basis ac- 



cording to 



8T 



e) = ^ ^ aimYern{0,4'), 



(38) 



=2 m=-e 



where 9 and (j) are the angles defining the direction along 
which the vector e is pointing. Then, the angular two- 
point correlation function can be expressed in terms of 
the multipole moments Ci as 



(39) 



and, as a consequence, the two-point correlation function 
of the temperature fluctuations operator can be written 



as 



5T ST 
1^(61)1^(62) 



00 



= — > '(2^+l)QP^(ei-e2), 

(40) 

the quantity Pg denoting Legendre polynomials. 

In order to pursue our demonstration that the CMB 
anisotropies are entirely determined by the quantum fluc- 
tuations, let us now express the multipole moments in 
term of the cosmological perturbation power spectrum. 
Upon using Eqs. (35) and (38), one obtains 



(41) 

and, from this expression, it is easy to show that 



1 



An f dk 



2a?Mlei 25 



-j1[fc(^o-r?&s)]7'.(fc), (42) 



where jg is a spherical Bessel function and where we used 
Eq. (34) to show that 



1 



2n^ 



T=>1 -'- 



V^{k)S{k-p). (43) 



We see that Cg is given by an integral over wavenum- 
bers of the Mukhanov-Sasaki power spectrum times a 
quantity that can be viewed as a "transfer matrix jik = 
3i (% ~ 'yfes)]" which allows us to "translate" a three 
dimensional spatial frequency fc into a two-dimensional 
spatial frequency I on the celestial sphere. We empha- 
size again that the above result is valid on large scales 
only; otherwise the integral in Eq. (42) contains another 
transfer function r^(fc) which takes into account the sub- 
sequent evolution of the modes when they re-enter the 
Hubble radius after inflation. Since ^ is a conserved quan- 
tity, we have T(;(fc — >• 0) = 1. 

Finally, let us also notice that Eq. (41) implies that 
{^\aira\% = since (^'ICfej*) = 0. Of course, this also 
means that l^^\ft IT\^) = 0. 
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E. Inflationary Predictions 

We have just seen that, in order to calculate the CMB 
multipole moments, we need to evaluate the curvature 
perturbation power spectrum. In this section, we calcu- 
late this quantity for power law inflation. 

The first step consists in solving the equation of mo- 
tion (12). Upon using Eq. (11), one obtains the time 
dependence of the frequency of the parametric oscillator, 
which reads 



a;2(77,fc) = k^ - 



/?(/? + 1) 
if' 



(44) 



From this expression, one sees that there are two regimes 
depending on whether the first term is dominant or sub- 
dominant. The Hubble radius is given hy l-^ = 1/ H — 
arj/(l -\- j3) and the Fourier mode wavelength can be 
expressed in terms of the co-moving wavenumber as 
A = 2'Ka/k. The first terms dominates if \kri\ ^ 1 or, 
equivalently, X <^ i-^. In this case lo k and we expect 
the mode function to oscillate as it would in Minkowski 
spacetime since, at those scales, spacetime curvature is 
negligible for the mode evolution. On the contrary, if 
|fcry| ^ 1, or A one has a; ^ l/ry, so curvature 

dominates and one obtains one growing mode and one 
decaying mode. These arguments are confirmed when 
one studies the exact solution for the mode function /fc. 
It can be expressed in terms of Bessel functions Ju{z) 
as [101, 102] 

fk = {-k-qf^ [CkJ/3+i/2i~kr]) + DkJ-(i3+i/2){-kr])\ , 

(45) 

where Ck and Dk are two integration constants. In or- 
der to match the initial vacuum behavior (24), one must 
choose 



TT C 



-i7r/4-i7r(/3 + l/2)/2 



2\l k sin[7r(/3 + l/2)] 



(46) 

In particular, one notices that both coefficients Ck and 
Dk scale as \j\fk. 

Since we want to evaluate the power spectrum on large 
scales, it is sufficient to take the limit fc?7 in Eq. (45). 
Then, one is led to 



1 



stand 2a2M2 £1 

PI ^ 

1 



7reim2/2 



/(/3)A:2'3+4^ Ag/fc^s-i, (47) 



where Mpj = TOpj/a/Stt and the fimction /(/?) is defined 
by [90] 



1 



r(-,g-i/2) 

21+/9 



(48) 



where r(z) is the Euler integral of the first kind [101, 102]. 
This function is such that, for the de Sitter case /3 = —2, 



one has /(/? = —2) = 1. The scalar spectral index 
rig =2/3-1-5 and, for solutions close to the de Sitter 
solutions, one has rig ~ 1, i.e. we have an almost scale 
invariant power spectrum. As discussed before, the devi- 
ations from scale invariance, are related to the deviation 
from the de Sitter case /3 = —2. This conclusion is in fact 
valid for any slow-roll models. The amplitude deter- 
mines the level of the temperature fluctuations observed 
in the sky namely 8T jT - 10"^ 

Finally, let us evaluate the multipole moments explic- 
itly. Upon using Eq. (42) and the expression of the power 
spectrum established above, one arrives at 

^ ^ ^^/^r[(3-n3)/2]r[^+(n3^1)/2] i_„^ 
^ r[(4-n3)/2]r[£ + 2-(n3-l)/2] ^ ^''^ 25' 

(49) 

where we have deflned r^ss = '7o~'7fss- Since this equation 
has been derived for large scales, roughly speaking one 
can estimate it to be valid in the regime I ^ 20. For ~ 
1, the above expression implies that Ci cx l/[^(^ -f 1)]- 

Of course, in the real world, the argument goes the 
other way round. From measurements of the CMB 
anisotropics, we observe that, on large scales, Ci oc 
l/[£(^ -f 1)] and, therefore, we deduce that the corre- 
sponding power spectrum is close to scale invariance, 
i.e. rig ~ 1. Obviously, this also means that a spec- 
trum that is not very close to scale invariance is now 
ruled out (more precisely, the WMAP data indicate that 



1 



= 0.018 



+0.019 
-0.02 



[20-22]). As already emphasized, 
the great success of inflation is that it precisely leads to 
such a power spectrum. 

It should also be clear that the above discussion, al- 
though perfectly correct at the level of principles, is over- 
simplifled at the technical level. The multipole moments 
are in fact computed at any scale (i.e. for any value of 
^) by means of numerical calculations (since, in the most 
general case, they are solution of more involved differen- 
tial equations) [103]. Moreover, their shape is not only 
determined by the spectral index but is also affected by 
the other cosmological parameters. The constraints on 
the different inflationary models are then obtained by a 
Markov Chain exploration of the parameter space [104]. 
But these technical considerations do not affect the con- 
siderations presented in this paper. Once again, as far 
as physical principles are concerned, the discussion pre- 
sented is this section is accurate. 



III. THE COSMOLOGICAL MEASUREMENT 
PROBLEM 

A. Squeezed State 

In this section, we study in more detail the properties 
of the quantum state in which the cosmological perturba- 
tions are placed [29, 31, 34, 105]. As already mentioned 
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around Eq. (19), it is described by the wavefunction 

^k{v,Vk,Vk) = ( — - — I e l"^"' i {50} 

^ ^QMenk^^ ^-2nMv^vl ^ (51) 

We see that this quantum state is completely known once 
the time dependence of rife (77) has been determined. The 
differential equation controlling the evolution of fl^ {rj) is 
given by the second of Eqs. (20). This equation is a Ri- 



catti equation (i.e. a first order, non-linear, differential 
equation). As is well known, it can always be reduced 
to a second order but linear differential equation. As 
already mentioned, this is achieved through the change 
of variable fl^ = ^ i/fc/(2//e)- The function fkii]) obeys 
fj! + uj^fk = and has been solved in Eq. (45). In the 
small-scale limit, one has fi^ — fc/2 and the wavefunc- 
tion (50) is the ground state of an harmonic oscillator. In 
the large-scale limit, a lengthy but straightforward cal- 
culation leads to 



k 



I 



(-/cr?) 



'r^-13-1/2) 



i-kri) 



-213-2 



(52) 



From this expression, one deduces that 



0, (53) 



and 



3mr2fc(7y) 



1 

2^ 



(l + /3) + ' 



2a 



— ^ 00, 



(54) 



where the limits are taken in the super-Hubble regime in 
which fcry — 0. 

We have mentioned above that the Ricatti equa- 
tion (20) can always be reduced to a linear second order 
differential equation. Of course, it can also be expressed 
as two linear, first order, differential equations. There- 
fore, one can introduce the functions Uk{ri) and Wfc('7) 
such that fk = (life -l- v^.) /V2k, the normalization 1/V2k 
being introduced for convenience. Then it is easy to show 
that these two functions obey 



Mfc — IKUk H . V 



(«y^)' 



(55) 
(56) 



The Wronskian W = f'kfl — Ik fk can be straightfor- 
wardly evaluated as W — i (|wfcp — |ffcP)- This means 
that, if we want to work with the choice W = i, one 
must have |ufeP — |wfcP = 1- This suggests to introduce 
the following parametrization 



Ukiv) = e"''=coshrfc, 



Vkiv) 



-'<^>^+^"t"' sinh 



rk- 



(57) 
(58) 



The three functions rk{ri), Okiv) ^i^d (f>kiv) ^^'^ called the 
squeezing parameter, rotation angle and squeezing angle 
respectively. It is clear that the knowledge of these three 



functions is equivalent to that of the function ^kiv) and, 
therefore, of the wavefunction. Upon using Eqs. (57) 
and (58), it is easy to show that 



cos {24>k) , (59) 
= /c - coth (2rfc) sin (20fe) , (60) 



(a /el)' 

k—- ^d— tanhrfc sin (20fc) . (61) 



The explicit relation between flk and the three squeezing 
parameters is given by 



k cosh Tfc — e sinh r'k a' 

i- 



2 cosh rk + e sinh rk 2a ' 
from which one deduces that 

k 1 

^e^k - 
5m flk = 



(62) 



(63) 



2 cosh (2rfc) -|- cos (20^) sinh {2rk) ' 
k sin (2^fc) sinh (2rfc) a' 

2 cosh (2rfc) -I- cos (2<^fe) sinh (2rfe) 2a' 

(64) 



Eqs. (59), (60) and (61) are highly non linear differ- 
ential equations and cannot be solved in general. We 
notice that Eqs. (59) and (60) are in fact decoupled from 
Eq. (61). Therefore, they can be solved in a first step and 
then the solutions can be inserted in Eq. (61) to find the 
behavior of Ok- In the case of power-law inflation, one 
can find explicit solutions for the de Sitter case, (3 = —2. 
Although this is not a solution for an arbitrary value of 
/3, it is sufficient to understand the main features of the 
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phenomenon of squeezing. One obtains 



rkifl) = -argsinh 



(t'hiv) = — I — arctan 1 , 
4 2 \2kf]J 



(2fcr;) ' 



(65) 
(66) 



Therefore, we see that, initiaUy in the sub-Hubble limit, 
rfc = (and 0^ = '^/'^) while the super-Hubble limit 
corresponds to the limit of strong squeezing — >■ +oo 
(and (f>k ^ 0). 

Based on the previous considerations, it is clear that 
the super Hubble limit is always associated with strong 
squeezing, even if we do not deal with the exact de Sitter 
solution. Indeed, now for an arbitrary /?, Eq. (60) can 
be written as (/)^ ~ — (/3 + 1) sin(20fc)/7y which can be 
integrated and leads to — arctan [C|77|^^(^+^)] . For 
/3 < —2, this confirms the fact that (pk 0. In the 
same limit, one has ~ l/rj from which one obtains 
Tfc cx (1 + /3)lna. This confirms that the super-Hubble 
limit is the strong squeezing limit and, given the fact that 
modes of astrophysical interest today leave the Hubble 
scale 50-60 e-folds before the end of inflation, one can 
deduce that — 120 for those modes [29, 30]. Compared 
to what can be achieved in the laboratory in quantum 
optics, this is a very large value [lOG]. 

In order to understand better the features of the quan- 
tum state (50), it is also interesting to calculate the mean 
values and dispersion of various quantities. First of all, 
it is clear that 



~R.I 



~R.I 
Pk 



0. 



Secondly, we also have 

2 



'k 



~R,I 
Pk 



(Smllfc) 



Finally, the cross-products can be expressed as 



\PkVk\ *) 



25Ref7fc' 

iflk 



-I + 



25Rel7fc' 



(67) 

(68) 
(69) 

(70) 
(71) 



and, of course, similar expressions for the operators 
v\ and p^. It is also interesting to notice that 

(*l«MI*) = (*l«l*> = o. ^ 

At this point, it is worth digressing about the 
definition of the conjugate momentum. The ac- 
tion (6) is of course defined up to a total deriva- 
tive. In Rcf. [15], it was shown that adding the term 
d {{d la)(v\f + {o! lo){v\f \ /(2d77) can also be viewed 
as a canonical transformation. This generates an addi- 



tional term (a'/°')(Pk '^^k 



R,I* R,I 



) in the Hamil- 



tonian. A complete study was presented in Ref. [15] 



and, here, we only quote the main results. It was shown 
that, at the quantum level, this canonical transformation 



R I 

leaves the amplitude ' invariant but induces the fol- 

-^R 

lowing transformations for the momentum: p^. 
with 



R,I 



.R,I 



-RJ 



-R.I 

Pk 



~R,I 



(72) 



On the other hand, the wavefunction is also modified, 
vl'fc — > ^fc, and the function £7^ changes according to 
Ofc — )• rjfc, where 



^k = ^k+i 



2a' 



(73) 



In particular, we see that the canonical transformation 
is such that the term ia' /(2a) in the expression (62) of 
the function ilkiv) is exactly canceled. The factor A'';, 
of the wavefunction is not modified and is still given by 
the first of Eqs. (21) (but of course should be used either 
with rifc or Hfc according to which set of variables is used) . 
This also means that when the averages (67), (68), (69), 
(70) and (71) are computed in the state \^), one obtains 
exactly the same expression, il^ being just replaced with 
Clk (of course, and \^), being related by a canonical 
transformation, represent the same physical state). 

We now come back to our calculation of the disper- 
sion of amplitude operator and its conjugate momentum. 
Upon using Eqs. (68) and (63), one obtains 



1 

2k 



[cosh (2rfc) + cos {2(j>k) sinh (2rfe)] 

(74) 

In the same manner, the dispersion of the operator tt^. ' 
is given by 



* = 



1 + sin2(2(/«fc)sinh^(2rfc) 



2 cosh (2rfc) + cos (20^) sinh (2rfc) 

(75) 

Let us now consider two new operators and Bp, ' , 

defined from tt^'Vv^ and Vkv^'^ through a rotation by 
the squeezing angle (fik- 

A^'^ = cos (j)k + ^Vk ^^ sin (/)fc, (76) 



-R,I 



Vk 



sm 



0fc - V kv^'^ cos 4ik- 



(77) 



It is easy to check that Ak,Bk = i- Then, a lengthy 
but straightforward calculation leads to 



A 



R,I 



B 



R,I 



(78) 
(79) 



Therefore, we see that there exists a direction in the plane 
(T^k^Vk) where the dispersion is extremely small. This is 
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why the corresponding state is called a squeezed state. 
In order to satisfy the Heisenberg inequality, the disper- 
sion along the direction perpendicular to the previous 
one becomes very large. As already mentioned, the phe- 
nomenon of squeezing is widely studied in many differ- 
ent branches of physics, in particular in quantum optics. 
Squeezing occurs each time the quantization of a para- 
metric oscillator is carried out. It is remarkable that the 
quantization of small fluctuations on top of an expand- 
ing universe also leads to that concept (squeezing here, 
i.e. Tfc ^ 0, does not require an accelerated expansion, 
only a dynamical background is necessary). 



ter Tfc is large. In the literature, this regime is very often 
described as a regime where the cosmological perturba- 
tions have classicalized [31, 32, 39, 39, 107]. Since this 
concept is subtle in quantum mechanics (and particu- 
larly when quantum mechanics is applied to cosmology) , 
we need to come back to this issue and to describe accu- 
rately what is meant by a "classical limit" in this context. 
In particular, it may seem strange at first sight that a 
quantum system placed in a strongly squeezed state can 
be described as a classical state since, in the context of, 
say, quantum optics, a similar situation would precisely 
be described as a non-classical situation [108, 109]. 



B. The Classical Limit 



We have seen in the last section that the super-Hubble 
limit corresponds to a limit where the squeezing parame- 



A convenient tool to study this question is the Wigner 
function, defined by 



X J 



■^,vi~ - ]e 



y 



(80) 



W{v^,vlp^,pi) = ^ I dxdy^* {vl 

Indeed, it is well known that the Wigner function can be understood as a classical probability distribution function 
whenever it is positive definite. Then, upon using the quantum state (50), the following explicit form is obtained 



27r3?e 



exp 



'23?er2t 



(p«-|-25mr!fc< 



exp 



1 

meOk 



(pj, + 23mr!fe4) 



(81) 



r 



The following remark is in order at this stage. One could 
have calculated the Wigner function with the state ^k- 
Obviously, one would have obtained exactly the same ex- 
pression except that all the fl^ terms would have been 

— Ft I • R, I 

replaced with fl^ and P)^' with tTj^' . In particular, this 
means that the term in parenthesis in the argument of 
the exponentials would have read tt^'^ -I- 25mr2fci;^'\ 
But, thanks to Eqs. (72) and (73), this is precisely 

R, I R I 

Pj^.' -I- 29mr2fcWj^' since the two terms proportional to 
a' /a exactly cancel out. This is of course related to the 
fact that the Wigner function is invariant under a canon- 
ical transformation. 

The Wigner function (81) is represented in Fig. 1 at 
different times or, equivalently, at different values of 
(rfc = 0.0005, 0.48, 0.88 and 2.31). The effect of the 
strong squeezing is clearly visible. Initially, in the sub- 
Hubble regime, is small and the Wigner function is 
peaked over of small region in phase space. As inflation 
proceeds, the modes become super Hubble and rf. in- 
creases. As a consequence, the Wigner function spreads 
and acquires a cigar shape typical of squeezed states. In 
fact, in the strong squeezing hmit, one has 3?ef^fc — >■ 
and 9mJ]fc — >■ fc sin 0^/(2 cos (/)fc) 0, see Eq. (63) 
and (64). Let us notice in passing that this last equa- 
tion is consistent with Eq. (54). On the other hand, if 



one considers 5m fife, then the leading term a' /{2a) is 
absent and one has to go to the next order in Eq. (64). 
This one is given by fc/[4(/3 -|- 3/2)]{—kf]) and represents 
the leading term of 5m Hj-. It goes to zero in agreement 
with the fact that 0^ ^ in the strong squeezing limit. 
In this regime, the Wigner function can be written as 



V'' COS^fc 



(82) 



This last equation represents the mathematical formula- 
tion of the cigar shape mentioned above. 

It is important to notice that the behavior described 
above is very different from the behavior of the Wigner 
function of a coherent state. The coherent states are 
usually considered as the "most classical" states and their 
Wigner function is given by 



(83) 



where v^''^^ and Pfc''^' represent the classical solutions. 
The typical shape is plotted in Fig. 2. One sees that 
the Wigner functions remain peaked over a small region 
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FIG. 1: Wigner function of a squeezed quantum state at different times during inflation. Only tfie two-dimensional function 
corresponding to the set of variables [v^,p^) has been represented, see Eq. (81). The time evolution of Kefifc and Sjmilfe 
has been expressed in terms of the two squeezing parameters rk and 4>k- These ones are given by the solutions (65) and (66). 
The left upper panel corresponds to rjj — 0.0005 and the corresponding state is almost a coherent one. The right upper panel 
corresponds to Vk = 0.48, the left bottom one to rk = 0.88 and, finally, the right bottom one to rk = 2.31. The effect of the 
squeezing and the cigar shape of Eq. (82) are clearly visible. 



in phase space and that this packet follows the classical 
trajectory (an ellipse in this context). Comparing Figs. 1 
and 2, we understand why a coherent state is usually con- 
sidered as classical while a squeezed state is considered 
as highly non classical. In the case of the coherent state, 
if one is given, say, the value of v^, then one obtains a 
value for the momentum, p^, which is very close to the 
one one would have inferred in the classical case. This is 
of course due to the fact that the Wigner function follows 
the classical trajectory and has minimal spread around 
it in all phase space directions. On the contrary, in the 
case of the squeezed state, if one is given then the 
value of is very uncertain since the Wigner function 
is spread over a large region in phase space. Therefore, 
we conclude that the cosmological perturbations do not 
behave classically in the usual sense. 

Given the previous discussion, it may seem relatively 
easy to observe genuine quantum effects in the CMB. 
Unfortunately this is not so, essentially because, in the 
strong squeezed limit, all quantum predictions can be 
in fact obtained from averages performed by mean of a 



classical stochastic process. 

Let us first study how this question is usually treated. 
For this purpose, let us consider again the expectation 
of the operator (tt^) [of course, one could also treat 
the case of ('^'fc )"]. The quantum average is given by 
Eq. (69), namely 



(84) 



(85) 



On the other hand, if one computes the quantity 



where W^rfc-j-oo (w^ , 'I'fc ) refers to the Wigner function 
in the strong squeezing limit (82), then one obtains 

(3mj7fe) /dieCtk, which coincides with Eq. (84) in the 
limit Tfc — >■ cx). This result is often taken as a proof 
that a strongly squeezed state can be described as clas- 
sical stochastic process. However, this argument is not 
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FIG. 2: Wigner function of a coherent state (83), represented at different times during inflation. Contrary to the Wigner 
function of a squeezed state of Fig. 1, the shape remains unchanged during the cosmological evolution. The Wigner function 
just follows the classical trajectory, an ellipse here since we deal with an harmonic oscillator. This justifies the fact that a 
coherent state can be viewed as the "most classical quantum state" . 



very convincing since it is a theorem [31] that the exact 
Wigner function [we stress again that, in Eq. (85), we 
have not used the general Wigner function but its limit 
when Tf, is large] satisfies the following property 



(86) 

where A is an arbitrary operator. Therefore, it does not 
come as a surprise that an expression like Eq. (85) repro- 
duces the corresponding quantum average in the limit 
Tk oo. 

In fact, as was discussed in Refs. [-34, 39, 110], what 
makes the situation so peculiar is something different. 
The point is that, in the limit — oo, all the quantum 
predictions can be reproduced if one assumes that the 
system always followed classical laws but had random 
initial conditions with a given probably density function. 
This can be easily understood on the example of a free 
particle [34, 39, 110]. Let us assume that, initially (at 
t = 0), the probability to find the particle at x is given 



by 




(87) 



where & is a parameter that characterizes the width of 
the distribution. At time t, this probability is given by 



1 



7r52 ^l + 4iV(m264) 



X exp 



26^ [x - kot/m) 
64 +4<2/m2 



where m is the mass of the particle and fcg the center of 
the Gaussian wave packet in Fourier space. 

Now let us consider a situation where we repeat many 
times an experiment consisting in sending a classical par- 
ticle from the origin with a velocity v (equivalently, in- 
stead of repeating the experiments many times, one could 
also consider an ensemble of classical particles) and de- 
tecting it at a position x ^ 0. By definition, the particle 
follows the laws of classical physics which means that its 
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motion can be described by the equation: x = vt (they 
all start from a; = at t = 0). Then, let us assume 
that the velocities are classical random variables with a 
probability distribution function given by 



P{v) 



1 



^/^^Av 



(89) 



This means, that according to the particle considered, 
the velocity is in fact not always the same. But because 
different particles have different velocities, they will not 
reach the position x at the same time. It is important 
to stress that, here, only the initial conditions are ran- 
dom and that the trajectory is purely classical. From the 
above distribution, we can easily infer that the probabil- 
ity of finding a particle at x, at time t, is 



Pix,t) 



1 



/ntAv 



-{x-vty /(tAvY 



(90) 



This distribution is in fact exactly |4'(x,i)p in the limit 
t — > oo provided we identify v — ko/ni and Av = 
\/2/ (mb) . Let us notice that this last relation is ex- 
actly what is obtained at the quantum level since x and 
V are conjugated variables. As a matter of fact, Eqs. (87) 
and (89) are Fourier transforms of each other. We con- 
clude that, provided we detect the particles far from the 
origin, the quantum predictions for the particles can be 
completely mimicked by means of a classical stochastic 
process. 

As discussed in Ref. [llD], the situation is exactly sim- 
ilar for the inflationary perturbations. The limit rp, oo 
is in fact equivalent to the limit of large times in the ex- 
ample above. One can even calculate the Wigner function 
of the free particle described by the wavefunction (88) 
and show that it takes the same form as the one of 
Eq. (81). Therefore, the inflationary perturbations are 
said to be classical in the sense explained before: they can 
be described by a classical stochastic process. In prac- 
tice, for instance, one can consider the aim in Eqs. (38) 
and (39) as classical random variables with probability 
density functions given by 



P{4ra) = 



1 p-«)V(2C,) 



2TrCe 



1 e-("-r/^', m^O. 



(91) 
(92) 
(93) 



Of course one can check that {aimai'm') = CiSwSmm' 
where, now, the bracket means a classical average calcu- 
lated by means of the above distributions. 

Finally, we conclude this section by a few words on the 
density matrix p^. In fact, the density matrix is nothing 
but the Fourier transform of the Wigner function. Let us 
denote by \v^') the eigenstates of the operator v^. Then, 



we have 

(94) 

Upon using Eq. (81) in the above equation, one arrives 
at 



xe 



(95) 



We notice that the off-diagonal terms, v^' ^ v^, oscillate 
very rapidly in the strong squeezing limit. This means 
that decoherence (defined as the disappearance of those 
off-diagonal terms) does not occur without taking into ac- 
count an environment for the perturbations. Various dis- 
cussions on what this environment may be can be found 
in Refs. [41-43]. 



C. Ergodicity 

Let us now discuss how, in practice, we can check the 
predictions of the theory previously reviewed. Initially, 
the system is placed in an eigenstate of the Hamiltonian 
(the vacuum state) of Eq. (19), which can also be ex- 
pressed as a superposition in the basis of the states \v^), 



namely 



1*)= / dz;« A^fc(,7)e-"'=('')K)V^) 



(96) 



The corresponding mean value of the Hamiltonian oper- 
ator can be expressed as 



1 



2 43?cf7fc 



(97) 



Of course, initially Ofc = k/2 and the energy is nothing 
but w/2 as expected for the vacuum state. 

In the real world, we measure the temperature 
anisotropics. As we have seen (and as is appropriate for 
an observable in the quantum- mechanical framework), 
this quantity is represented by an operator. According 
to Eq. (38), measuring the temperature anisotropics is 
equivalent to measuring the observables which, in 
turn, according to Eq. (41), is equivalent to measuring 
the observables Cfc or (that is to say and v\). 

According to the postulates of quantum mechanics, 
measuring the observable ii^ gives an eigenvalue (no 
hat, it is a number) with probability K'l'fe l^*)!^ and, im- 
mediately after this measurement, the system is placed in 
the eigenstate More concretely, after the measure- 

ment, we "see" a specific CMB map and we say that the 
measurement has produced a specific "realization". The 
result is given in terms of coefficients ai„i (again, no hat) 
expressed in terms of the numbers through Eq. (41) 
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(except, of course, that this equation should now be used 
with no hat on both sides). Equivalently, we see a specific 
temperature pattern ST{e)/T (no hat) corresponding to 
the set of numbers aim , see also Eq. (38) . In conclusion, 
the CMB map observed, say, by the WMAP satellite cor- 
responds to one measurement (or one "realization") of 
the operator dT{e)/T. 

Then comes the question of how one can operationally 
verify these theoretical predictions. In quantum mechan- 
ics, in an ordinary laboratory situation, one would check 
that the theory is correct by repeating the experiment 
many times. In this way, one would generate many real- 
izations of (or, equivalently, of a^m. or ST/T) i.e. one 
would obtain iVVcai numbers uj^^ z = 1, • • • , A^rcai [or 
or {ST/TY] where A^rcai is the number of realizations 
(that is to say the number of times the experiments have 
been performed). With these A^roai CMB maps, one could 
then check that the v^'' are indeed distributed with a 
Gaussian probability density function in agreement with 
Eq. (50) or, with the A^^^eai sets of numbers a^^, one could 
infer whether they follow Eqs. (91), (92) and (93), de- 
termine the corresponding variance and check that it is 
given by the Ce predicted by the theory. Let us notice 
that the above discussion is independent from the fact 
that the perturbations can be described classically or not. 
If we are in the classical limit (in the restricted sense de- 
fined in the previous section), then we showed that mea- 
suring the observable can be viewed as measuring a 
classical system with random initial conditions but this 
does not change the fact that we need many realizations 
to check that the probability density function predicted 
by the theory is the correct one. 

Clearly, in cosmology, the program described above 
cannot be carried out because one cannot repeat the ex- 
periment many times since we are given only one CMB 
map [-V^]. How, then, can we check the predictions of 
the theory of cosmological perturbations? To discuss 
this question, let us be more accurate about the oper- 
ator ST/T{e). In the large-squeezing limit, we have seen 
that it can be viewed as a classical stochastic process and, 
therefore, it is convenient to write it as 

f^(C,e), (98) 

where the symbol ^ labels the realizations. A given re- 
alization of a stochastic process is a function of e. By 
contrast, a given realization of a random variable is not 
a function but a number. This is for instance the case 
of a^m(^). The idea is then to replace ensemble averages 
by spatial averages (i.e. averages over different directions 
e) [33]. If the process is ergodic, these two types of aver- 
ages are equal [33]. In that case, one can check the pre- 
dictions of the theory even if one has only one realization 
at our disposal. Unfortunately, one can also show that 
a stochastic process living on a sphere (here, of course, 
the celestial sphere) cannot be ergodic [33]. Therefore, 
we are left with the task of constructing unbiased estima- 
tors with minimal variances. For instance, let us assume 



that we have calculated the number Ce in some inflation- 
ary scenario and that we would like to compare its value 
to an actual measurement. How would we proceed? We 
would consider the random variable Ce{(,) defined by the 
following expression [33] 

CiiO^i- I df]idf^2P^(cos<5i2)^(e,ei)^(^,e2), 
47r T T 

(99) 

where 5i2 is the angle between the direction ei and 62. 
As announced, the estimator C^(^) is expressed as a spa- 
tial average of the stochastic process 5T/T. It is easy to 
show that it is unbiased, {{Ci)) = Cf and has the min- 
imum variance [33] (called the "cosmic variance") given 
by y/2/{2l^\)Ci>. The double brackets (()) means an 
ensemble average, which amounts to a quantum average 
in the high squeezing limit as mentioned before. One 
should be careful that this ensemble average has noth- 
ing to do with the one introduced below (denoted E) 
for the CSL modifications of the Schrodinger equation, 
since these two stochasticities have completely different 
natures, the former being effective and the later intrinsic. 

In practice, we would proceed as follows. From 
our CMB map ST{^,e)/T, we compute the integral in 
Eq. (99) and this gives a number representing one realiza- 
tion of the estimator Ci, the only one we can have access 
to. It is unlikely that this number be because it is un- 
likely that one realization of a random variable be exactly 
equal to the mean value of that variable. However, if the 
variance is small (i.e. if the estimator is good), the cor- 
responding probability density function will be sharply 
peaked around the mean value and any realization will 
therefore be close to the mean (and, in our case, it is 
not possible to decrease the value of the variance since 
we work with the best estimator). Therefore, we can 
study where the number we have obtained by following 
the above described procedure falls, compared to the in- 
terval Ci ± y/2/{2£ + l)Ce, where Ce is the theoretically 
predicted multipole moment. Then, for instance, one can 
start a calculation of the to assess to which confidence 
we have verified the theory. In fact, the cosmic variance 
can simply be seen as another source of error, besides 
those coming from the instruments. 

Given the previous discussion, there is one issue that 
one can raise and which is the subject of the present 
paper. The question is how a specific outcome (a real- 
ization) is produced. Above, we have just assumed that 
this happens without discussing this point. According 
to the postulates of quantum mechanics in the Copen- 
hagen interpretation, this "macro-objectification" takes 
place when a measurement is performed. Since the CMB 
anisotropics were produced at a redshift of ziss — 1100, 
this means that it should have happened prior to that 
epoch (possibly during inflation itself). But, clearly, 
there was no observer at those early times. We face 
here the conventional measurement problem of quantum 
mechanics which is, in the context of cosmology, exacer- 
bated. 
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IV. THE PEARLE-GHIRARDI-RIMINI- WEBER 
THEORY 

A. A Dynamical Collapse Model 

Although one can manage to obtain, based on primor- 
dial vacuum quantum fluctuations, a set of correlation 
functions that are formally indistinguishable from a clas- 
sical stochastic distribution, one still has to face the prob- 
lem of reaching a specific realization before cosmological 
perturbations can start to grow in a classical way. This 
amounts to the question of the measurement problem in 
quantum mechanics, namely that there are two distinct 
evolution processes: the unitary and linear Schrbdinger 
time evolution on the one hand, and the stochastic and 
non linear wavepacket reduction on the other 

In what follows, we briefly present the collapse the- 
ories and explain how the Schrbdinger equation can be 
modified in order to allow a dynamical description of the 
wave-packet reduction. In fact, to be more precise, we 
shall restrict attention to the case of CSL [G2, 63, 65, 67]. 

The CSL model relies on the idea that to the 
Schrodinger linear evolution should be added an extra 
stochastic behavior, encoded through a Wiener process 
whose differential acts as a random square root of 
that of time, namely 

E{dWt)^0, and E {dWtdWf) = 6 (t - t') dt^ , 

(100) 

where E stands for an ensemble average. One then ex- 
pands the state vector variation d|x) up to first order in 
time through 



d|x) = (^Adt + BdWt) \x), 



(101) 



where A and B are operators acting on the Hilbert space 
of available states. One then demands that, on average, 
the wavefunction be normalized, i.e. 



E((xlx)) = l 



E[d((x|x))] =0, (102) 



which, upon using Ito calculus^ . for the differentials and 
Eq. (100), yields 



+ i = -B^S, 



(103) 



since the state |x) is arbitrary. The general solution of 
Eq. (103) is A = —iH — ^B^ B, where H is hermitian 
and to be identified with the Hamiltonian leading to the 
usual Schrodinger dynamics. 



^ This means that for two functions / and g of the stochastic 
variable W, one has d(/g) = fdg + (df)g + E [(d/) (dg)] and 
df{W) = f'{W)dW+y"{W)E [(diy)^] , where a prime stands 
for ordinary derivative with respect to the argument W. It is 
necessary to expand up to second order in the noise because 
Eq. (100) means E{dW^) 



In order to assign a probabilistic meaning to the norm 
of the wavefunction, it should be normalized. However, 
according to Eq. (101), although this is true on average, 
it varies stochastically according to 

dllxll' = (Xl {b + B^) \x}dWt - 2{x\B\x}dWt, (104) 

where from now on we assume that B is hermitian. 

Eq. (104) implies that the state |x) is not normalized, 
and one can define a normalized one whose probability 
distribution will thus be interpretable in terms of mea- 
surements. We then set 



11x11' 



(105) 



whose dynamics can be computed using the previously 
derived rules. One finds 



d|V) = 



-^H--[B-{B] 



+ (B-{B))dWt\\i;), 



dt 



(106) 



where the quantum expectation value is taken on the 
normalized state vector and thus defined as 



{B) = {^\B\^). 



(107) 



The operator B can be decomposed as B — ^Q. The 
coupling constant 7 is the product of the localization rate 
with the width of the Gaussian wavefunction inducing 
the localizations [62], and sets the strength of the non- 
linear effects and therefore the characteristic time scale 
over which these are measurable. The observable Q, for 
instance the position operator, is the basis on which the 
states are to spontaneously collapse to (in the following, 
we also call the operator Q, the "collapse operator"). 

As it turns out, and this is exemplified later in the case 
where the operator Q is identified with a cosmological 
perturbation Fourier mode (see Sec. VA), the natural 
evolution of Eq. (106) is to project an initial state I^/iq) 
on an eigenstate |a) of the operator Q: setting 



■■ dt. 



(the sum being replaced by an integral in the case of a 
continuous spectrum for Q) such that Q\a) = qa\a), one 
finds that limt_>.oo |^(i)) — \ct) for a given value of a, and 
this with a probability P{a) = |(^'|a)p. In other words, 
the Born rule is naturally implemented as a dynamical 
consequence instead of being imposed as an extra hy- 
pothesis. 

Finally, defining the density operator as 

p = Em{^\), (108) 

one obtains, using Eq. (106) the so-called Lindblad equa- 
tion, namely 

d/5 ■ r/y -1 7 



Q,P 



(109) 
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providing its time development. 

Let us now come to another very important aspect of 
the CSL theory and describe the so-called "amplification 
mechanism" which enables to understand why the dy- 
namics of microscopic systems is not much altered by the 
extra stochastic and non linear terms in Eq. (106). This 
is phenomenologically very important since this means 
that the laboratory experiments performed on "small" 
quantum systems are still accurately predicted by the 
standard Schrodinger equation while the macroscopic ob- 
jects are quickly and efficiently localized. Let us consider 
an ensemble of N identical particles, assuming that, for 
each of them, the collapse operator is the physical po- 
sition in space. Therefore, we can identify the operator 
and Wiener processes according to 



N 



B 



and dWt dW. 



(110) 



in Eq. (106), with Xi the position operator for the i*'^ 
particle. Note that in this case, one has as many in- 
dependent Wiener processes as there are particles; they 



satisfy 



E 



dwy^dwjp 



5'^ {t~t') df 



(111) 



This naturally generalizes Eq. (106) to a set of operators 
and particles on which to project the relevant states. 

We now assume that one can decompose the total wave 
vector l^f) in the form 



I*({^J)) = I«'CM (i?))^|*rel({rj)) 



(112) 



where the total wavefunction depends on the set of all the 
position operators {xi\, while the "macroscopic" part of 
it, I^'cm); depends only on the position R = N^^ Xi 
of the center of mass, and the rest is a function only of 
the relative coordinates defined through Xi = R + r^. 

Using Ito calculus to evaluate the differential of the 
tensor product in Eq. (112), it is easily checked that 
|*({a;J)) satisfies Eq. (106) with B and dWt given by 
Eq. (110) if the components of the product respectively 
satisfy 



and 



d|*cM {R)) = 



dt + V7^(i?-(i?))dW^t||vE',„(i?)), (113) 



T \ ^ 



i=l 



N-1 



+ V7 E - (^^)) dW^t {{n})) , (114) 



1=1 



where we have assumed the total Hamiltonian could be 
split into H = H^^^ {R)+H,-ci {{fi}) and the new constant 
7cM appearing in Eq. (113) is given by 7(,„ = N-y. This 
illustrates the mechanism thanks to which localization 
is amplified for a macroscopic object containing a large 
number (in practice N ^ 10^^ ^ 1 for usual classical 
systems) of particles, while the usual quantum spread is 
mostly conserved for the internal degrees of freedom. A 
recent inventory of all the constraints derived so far in 
various physical situations on the CSL parameter 7 can 
be found in Ref. [111]. 



An Illustrative Example: 
Oscillator 



the Harmonic 



in the case of the CSL theory. Moreover, in cosmology, 
as explained before, we deal with a parametric oscillator, 
a case which shares some similarities with an harmonic 
oscillator, at least in some regimes. It is therefore im- 
portant to understand first this simplest case in the CSL 
framework. In the following, we assume that the opera- 
tor B introduced in the previous section is the position 
operator x. As a consequence, the modified Schrodinger 
equation can be written as 



d* = 



--Hdt + ^-f{x-{x))dWt 



l{x-{x)fdt 



(115) 



In this section, we illustrate how the CSL theory works 
on the example of the harmonic oscillator resetting the 
Planck constant H for easier comparison with previous 
works. This is an interesting case because it represents 
the prototypical example of a quantum system and, to 
our knowledge, this case has not been solved explicitly 



where H — p^/ {2m) + mixP'5?' j2 is the Hamiltonian. The 
parameter 7 sets the strength of the collapse mechanism 
and, since we have chosen the position as the preferred 
basis, it has dimension x T^^ . Following Ref. [7-')], 
the wavefunction can be taken as a Gaussian state and 



19 



the most general form can be expressed as 
'^(t,x) = \N{t)\exp^^^^en{t)[x-x(t)Y 



+ix{t)x - 



i^mn{t)x^ 



ia{t) 
(116) 



where, a priori, |iV|, Sfteil, x, a, x ^i^d 9m f2 are real 
stochastic variables. Introducing this wavefunction in 
Eq. (115), one obtains the following set of equations 



\N\ 
(3?ef7)' 

(3m 17)' 



t:;^ = — 3ni 

4 sRel7 m 

7 + 4— (SRefl) (3m f7) , 
m 



7 



(117) 
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2{^enY - 2(Q:mr2)^ 



TO 

7i 2 



(118) 

2 



2(5mf7)x] 



V7 dWt 
2^cVL dt 



-5Re f7 + 2 (sRe f7) - 



1 



(119) 
(120) 

(121) 
(122) 



where a prime means a derivative with respect to time. 
We see that the first equation can be integrated to give 
|iV| = (23?efl/7r)i/4, which ensures that the wavefunction 
is properly normalized. Then, the two following equa- 
tions, Eqs. (118) and (119) "decouple" from the other 
equations and can be integrated separately. In particu- 
lar, if we add them up, we arrive at 



n' 



m 



(123) 



This equation should be compared to Eq. (21). As ex- 
pected, there are identical provided we take h ~ m ^ I 
and 7 = 0. Of course, in the present case, the frequency 
to is constant since we deal with an harmonic oscillator 
rather than a parametric oscillator as it is the case for cos- 
mological perturbations. Eq. (123) is a Ricatti equation 
and we have already seen that the appropriate change of 
variable to transform it into a linear second order differ- 
ential equation is = —imf'/{2hf), where the function 
/(<) obeys the equation 



uj^ - 2i— 7 
m 



/ = o. 



(124) 



This equation admits simple solutions that can be 
expressed in terms of exponentials, namely f{t) oc 
exp (±a t) where a is defined by 



2ijh 



(125) 



As a consequence, the solution for n{t) can be written as 



n{t) = — — atanh (at + . 



(126) 



where (/) is an integration constant that can be expressed 
in terms of the initial value of the function D,{t) 



argtanh 



2hn{t^ 0) 



(127) 



This solution resembles the formula obtained in the case 
of the free particle, see Ref. [73]. 

At this stage, we need to discuss the initial conditions. 
Our assumption is that, at < = 0, the quantum state is 
simply given by the ground state of the harmonic oscil- 
lator in conventional quantum mechanics. Technically, 
this means that we require the wavefunction to be 



v,(t = o,x)=f^V^'e— ^A^'^) 



V irh J 



(128) 



which implies that Jfteil = muj/{2h) and 5m = or, 
equivalently, (j> = argtanh(ia;/a). Notice that this choice 
is fully compatible with the normalization established 
above, |A^| = (23?ef^/7r)i/4. of course, our choice also 
amounts to imposing x{t = 0) = a{t = 0) = x(^ = 0) = 
0. 

Since the evolution of the stochastic wavefunction is 
controlled by the function ^l{t), it is interesting to study 
how it evolves with time. Writing the number a as a = 
+ ia^, where it is easy to show that 



a 



a 



and (j) = 0^ + i't'^i straightforward algebraic manipula- 
tions lead to the following expressions for 3?e and 3m fl: 




TO sinh [2 {aH + (f>^)] + sin [2 {aH + (f>^)] 

2h cos [2 {aH + 0i)] + cosh [2 (a^t + (j)^)] 

TO sin [2 {ah + (j)^)] - sinh [2 {a^t + (j)^)] 

2h cos [2 {aH + 01)] + cosh [2 {aH + 0R)] 



(131) 
(132) 
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In particular, the function ^e^(t), with the initial condi- 
tion specified above, is always positive. Notice also that 
there is a sign ambiguity in the definitions of the quan- 
tities and in Eqs. (129) and (130) but one can 
show that this does not affect the physical predictions of 
the model. It is also interesting to calculate the limit for 
large times of the two functions in Eqs (131) and (132). 
One obtains 



lim Kerj 

t—>oo 



lim 5m D, 



2h 



ma 



2h 



1 



1 te-f^ 



2 m2w4 



2H 



1_ 



1 h^-i 



2 m'-uj^ 
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■)• 

(134) 



where the dots indicate an expansion in the small di- 
mensionless parameter hqlimuP'). We see that, if 7 = 0, 
we obtain the ground state given by Eq. (128). Devia- 
tions from that solution are controlled by the parameter 



hqj (tooj^). 

We are now in a position where one can investigate 
the physical properties of the quantum state (116). In 
particular, it is easy to show that (i) = x and (p) = 
X — 2 (3m ri) X. Initially, x — ^ and the position operator 
has a vanishing mean value as expected for the ground 
state of the harmonic oscillator but, at later times, due to 
the stochastic evolution of the wavefunction, it acquires a 
non zero value. It is also possible to calculate the spread 
in position and momentum. One obtains 



1 1 



2 Vii^' 



(135) 



(5Rel7)2-h(!3ml7)2 



3?ef7 
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We see that these quantities only depend on SRefi and 
3m f2. As a consequence, inserting Eq. (131) and (132) 
in the above expressions of Ox and Up, one arrives at 




cos [2 {aH <^^)\ -h cosh [2 {a^t -I- 0^)] 
I sinh [2 (o^t + sin [2 {aH + (/)i)] ' 



cosh [2 {a^t + 0R)] - cos [2 {aH + ^^)\ 
sinh [2 [a^t -t- 0^)] -f a^ sin [2 (aH + (/)i) 



(137) 



The time evolution of these quantities is displayed in 
Fig. 3. The black curves correspond to the conventional 
Schrodinger evolution, i.e. the case 7 = 0. They show the 
usual oscillatory behavior. On the contrary, when 7 7^ 0, 
we see that the oscillations are damped (see the red and 
blue curves). Then, the spreads converge towards a con- 
stant value, which only depends on 7, to and m. This 
value is easy to evaluate and one finds 



lim cTj. = ^777 . /— 




1/4 



1 + 4^-1) ,(138) 



/ 2 ^ 1/4 

lim cr„ = muj 1 + 4 — - — 7 I x lim a,-- (139) 
t-s-oo P \ m^uj'^ J ^ ' 

From these formulas one can see that the spread in posi- 
tion at infinity decreases with 7, from ^JKJJ^moJ) for 
7 = to for 7 — >■ 00. We see that the modified 
Schrodinger equation, as expected, implies a localiza- 
tion in position. We also notice that the microscopic 
behavior of the system is altered by the non linear and 
stochastic terms added to the theory. By contrast, in 
order to satisfy the Heisenberg uncertainty relation, the 
spread in momentum increases with 7, from ^ mLuh/2 
for 7 = to infinity for 7 — > 00. For 7 = and at large 
times, one finds that the Heisenberg relation is saturated. 



o'xO'p = ^/2, as appropriate for a coherent state. In the 
limit 7 — cx) , one finds a larger value UxCTp = 
Let us also remark that an exact eigenstate of the oper- 
ator X is given by a Dirac function 8(x — x) centered at 
some value x. On the other hand, we see that adding non 
linear and stochastic terms results in a spreading of the 
Dirac function into a Gaussian wavefunction with a finite 
width decreasing for increading 7. Therefore, the mod- 
ified Schrodinger equation does not exactly lead to an 
eigenstate of the position operator. In fact, the asymp- 
totic value of Ox obtained above defines the "precision" 
of the collapse and characterizes how close to an eigen- 
state of the collapse operator the final state is. In that 
sense, since Ox decreases with 7, the bigger 7, the more 
"precise" the collapse. 

To conclude this section, it is also interesting to cal- 
culate the time derivative of the quantum mean value of 
the Hamiltonian operator. One obtains 



di 



2m^ m 



^^mfl dWt 



SRefl 
V7 dm 
SfteO At 



(140) 
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FIG. 3: Spread in position and momentum for different values of 7 in the case of the harmonic oscillator, see Eq. (137). 
The conventional Schrodinger evolution corresponds to 7 = and is represented by the black curve which oscillates. On the 
contrary, when the collapse mechanism is turned on, the oscillations are damped (blue and red curves), the spreads tend toward 
a constant value and localization occurs. 



This equation implies that 
dE (^H 



2m 



7- 



(141) 



As is well know, this formula expresses the non conser- 
vation of energy in the CSL theory. From a phenomeno- 
logical point of view, this increase of energy is usually so 
small (given the values of 7 usually considered) that it 
cannot be detected. Put it differently, the non conserva- 
tion of energy in the CSL theory cannot be used to rule 
out this theory [OS]. 



V. THE INFLATIONARY CSL THEORY 

The dynamical collapse model of the previous sections 
should apply to any quantum system, and hence in par- 
ticular to cosmological perturbations as they arise from 
vacuum fluctuations. Spontaneously collapsing these 
happens to be a tremendously complicated task for many 
reasons discussed below, so in what follows, we suggest a 
much simplified modeling method which we then apply 
to the inflationary situation. 



The Modified Schrodinger Equation for the 
Mukhanov-Sasaki Variable 



The first obvious problem one encounters when deal- 
ing with quantum cosmological perturbations is that the 
underlying theory ought to be relativistic. The straight- 
forward relativistic generalization of the CSL model for 
quantum field theory, starting with the action (6) in the 



Tomonaga picture for instance, leads to unremovable di- 
vergences [(>•")] (see however [G9]), even more so when non 
linearities inherent to general relativity are taken into 
account. 

The second next option which happens to lead to a 
model in which calculations are actually possible con- 
sists in noting, as mentioned earlier in Sec. II D, that 
the spectrum of primordial perturbations depends on the 
wavenumber k. In other words, once the Fourier spec- 
trum is known, all the observable quantities related with 
the CMB can be computed and compared with actual 
data. This means that mere knowledge of the modes Vk 
ought to be enough in order for a complete description 
of the possible observations to be realized. 

We shall therefore accordingly assume in what follows 
that the modified Schrodinger equation of motion for the 
wave function will be done at the level of the Fourier 
mode 'ilk , with spontaneous localization on the Vk eigen- 
manifolds. This is consistent with previous approaches 
aimed at studying decoherence of cosmological pertur- 
bations where the pointer basis is often assumed to be 
precisely the Mukhanov-Sasaki operators, see Ref. [38]. 
Separating as before into real and imaginary parts, we 
shall thus assume the following basic equation 



k J 



(142) 



and a similar equation for 5*^. Here, the quantity 7 is 
a positive constant with mass dimension 2 if the scale 
factor is chosen to be dimensionless but is dimensionless 
if the scale factor is chosen to have mass dimension — 1, 
which is the convention adopted here. As in Sec. IV, the 
parameter 7 sets the strength of the collapse mechanism. 
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FIG. 4: Evolution of various physical length scales with time during the cosmic history in the CSL model described by Eq. (142) 
with a zoom on the transition from inflation to reheating inserted (see the concluding section). The solid line represents the 
Hubble radius l-^ and the dashed-dotted green and red lines, the physical wavelengths of two Fourier modes of cosmological 
relevance today. The solid blue line represents the built-in CSL scale i-y, see the discussion above Eq. (156). It is a preferred 
comoving scale and can also be viewed as a time dependent preferred physical scale. Therefore, when a mode is below (above) 
l~f it remains so during the whole history of the universe as is clear from the plot. This means that, contrary to the Hubble 
scale, there is no "^.^ crossing" during the cosmic evolution. As a consequence, one expects the power spectrum to acquire a 
broken power-law shape, with two different branches, an expectation confirmed by the calculations in Sec. VD. 



Let us now review all the limitations of postulating an 
ansatz equation such as Eq. (142). First, one should note 
that the constant 7 in Eq. (142) cannot be the same as the 
one associated with the choice of the position operator as 
the collapse operator appearing in Eq. (115), despite our 
choice of the same notation. It is clear that each time 
one considers different "collapse operators" , this leads to 
different CSL parameters with different mass dimension. 
The same phenomenon is observed in Ref. [70] where the 
"collapse operator" is chosen to be a spin operator. In 
this case, it is clear that the corresponding CSL param- 
eter cannot be the same as the one corresponding to the 
case where the "collapse operator" is the position (as it is 
for the case of the free particle [7.S]). This is unfortunate 
when it comes to a comparison of the constraints ob- 
tained in the laboratory with the constraints obtained in 
cosmology. In fact, what could be done is to consider the 
strict CSL theory where the "collapse operator" is usu- 
ally taken to be the averaged density operator. In the 
language of cosmological perturbations, this amounts to 



assuming that there is spontaneous localization on the 
(5/9(77, x) eigenmanifold, where 5p{ri, x) is the perturbed 
energy density. This would have the advantage to in- 
troduce a universal 7 with always the same dimension. 
Unfortunately, 5p{rj, x) is a complicated functional of 
and this would probably render the whole approach un- 
tractable. Let us also notice that 7 could be taken as a 
function of the wavenumber fc, i.e. different CSL param- 
eters for different modes. In this article, for simplicity, 
we do not follow this route. 

Another issue is that we moved from real to reciprocal 
space while keeping the structure of the equation un- 
changed. In so doing, we also avoid from the outset any 
mode mixing that would be naturally arising from a real 
space modified Schrodinger equation: its stochastic ver- 
sion being non linear, one would expect a coupling of the 
Fourier modes, which is here automatically set to zero. 
Note this approximation is justified by data observations 
of the CMB. 

Another important limitation of our treatment is the 
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fact that the collapse concerns the modes independently. 
As a result, the amplification mechanism, so crucial to 
explain why the quantum behavior becomes increasingly 
less important for increasingly large systems (the effec- 
tive collapse time being inversely proportional to the 
number of particles involved and, hence, to the size of 
the system), is simply not operating here! Therefore, 
even though one might consider cosmological size effects, 
the collapse will occur just as it would for an indepen- 
dent quantum particle. As we will see, that implies a 
severe constraint on the constant 7 when comparison of 
the modified spectrum is made to actual observations on 
Hubble-size scales. 

Finally, Eq. (142) is written in terms of the confor- 
mal Fourier mode of the original action. Because its 
normalization implies the equation be non linear, this 
means the constant 7 can be translated, as we will show 
later, into a privileged conformal scale, and hence a time- 



dependent privileged length i^, as shown in Fig. 4 and 
the discussion above Eq. (156). This is somehow similar 
to the fact, except at the perturbative and conformal lev- 
els, that considering non flat spatial sections permits to 
define a curvature length and thus forbids to renormalize 
the scale factor arbitrarily. However, as shown in Ap- 
pendix A, this last limitation does not affect the general 
conclusions drawn here. 



B. Gaussian State 

Our goal is now to solve Eq. (142). As was done for the 
standard case (50), one considers that the wavefunction 
assumes a Gaussian shape. Concretely, we take the most 
general form, namely 



RJ 



-R,I 



iv) 



1(7,. 



RJ 



R,I 



i3m fl 



Uv){4-'y}, (143) 



where v 



R,I 
k 



R.I 1 n 

cr,. and Xk 



R,I 



are real numbers. The fact that one can assume \Nk\ and to be independent of "R" 
or "I" will be justified in the following. Compared to Eq. (50), we see that Eq. (143) is more general and, therefore, 
contains more parameters. The case of Eq. (50) corresponds to v^'^ = 0, Xfc'^ = and arg A^^ = Cfc. Of course, the 
above Gaussian is similar to the wavefunction considered in the case of the harmonic oscillator of Eq. (116). The only 
difference is that the stochastic functions characterizing the wavefunction now depend on the wavenumber fc and the 
role of the position is played by the Fourier amplitude of the Mukhanov-Sasaki variable. 

The next step is to insert Eq. (143) into Eq. (142) in order to derive the differential equations obeyed by the 
functions parameterizing the Gaussian state. Straightforward manipulations making use of the Ito calculus lead to 
the following expressions 



\Nk\' HJRenkY ^ „ , 7 
|A^fc| A ^cQk 43ficOfc' 

(3?ef7fe)' = 7 + 4(5Rel]fe)(3?nif^fc), 

Qk)' = ~2i^enkf + 2{^mnkf + ^uj^iTi,k) 



(3m 



(' 

(Xfe 



R,I 



R„I 



= -3?e rtk + 2 (5Re ilkf {v^^') ' - ^ {xk'') ' , 
= -4(3fier!fe)'f;^'^ + 2xfc''(3mr!fc). 



(144) 

(145) 
(146) 

(147) 

(148) 

(149) 



Several remarks are in order at this point. Firstly, we see 
that the evolution equations for lA^^I, ^eQk and 3m fife 
are deterministic and independent of that of v^'^, a^'^ or 

R I 

Xk ■ This justifies the fact that one can assume these 
quantities to be independent on R, I provided similar ini- 
tial conditions are chosen for R, I. This also means that 
these three quantities are not random (but their evolu- 
tion is still explicitly modified by the stochastic dynamics 



when 7 7^ 0). Secondly, Eq. (144) explicitly implies the 
conservation of the wavefunction norm: if one initially 
has a normalized state, i.e. 

m^{^Y\ (150) 

it will remains so at any time. In fact, this equation is 
similar to Eq. (21) which is therefore not modified by the 
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introduction of the non hnear stochastic terms. More- 
over, in the present case where the wavefunction is given 
by a single Gaussian, a^'^ is just an irrelevant global 
phase and can be ignored (this will no longer be the case 
when the quantum state is the sum of two Gaussians, see 
below). Thirdly, it is easy to check that Eqs. (144), (145), 
(146), (147), (148) and (149) are the exact counter parts 
of Eqs. (117), (118), (119), (120), (121) and (122). The 
only difference is that uj is now a time-dependent quan- 
tity as expected since we deal with a parametric oscilla- 
tor. We conclude that, instead of six coupled stochastic 
differential equations, we have in fact to solve two sets 
of two coupled differential equations, the first one be- 
ing deterministic and the second one being stochastic. 
In particular Eqs. (145) and (146) can be combined and 
lead to the following Ricatti equation for the quantity 
3?e rifc -|- iQm fl^ = i^k - 



(151) 



This equation is similar to Eq. (123) obtained for the 
harmonic oscillator. Of course, if 7 = 0, then one ex- 
actly recovers the Ricatti equation (20). As discussed be- 
fore, a Ricatti equation can always be reduced to a linear 
but second order differential equation: this is achieved 
through the transformation flk = —ifk/C^fk), where fk 
is a solution of the following linear differential equation: 



f'^+[Lu' (77,fc)-2*7]/fc=0. 



(152) 



This equation is very similar to the equation for the mode 
function considered before. The only difference is the 
appearance of the term —2^7 in the effective frequency. 
Obviously, if 7 = 0, then one recovers the conventional 
case. Moreover, the fact that this is the counterpart of 
Eq. (124) is obvious. 



C. Evolution of the Stochastic Wave- function 
during Inflation 

We now study the time evolution of the quantum 
state (143) in more detail. We start with the evolu- 
tion of 5fte Qk and ^Ik since we have shown in the 
last section that it decouples from the other equations 
of motion. To derive the corresponding solutions, it is 
sufficient to solve Eq. (152). If the background is driven 
by a phase of power-law inflation, w (77, k) is given by 
UJ (77, k) = k"^ — P [P + \) / rf' and the differential equa- 
tion (152) reads 



fk 



/3(/3 + l) 



77^ 



- 2i7 
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We see that the only effect of the CSL term —217 is to 
modify the comoving wave number fc^ — > fc^ — 217. The 
solution of Eq. (153) can be written in terms of Bessel 
functions 



fkivi) = i-Zkkriy^^ [CkJp+i {-Zkkrj) 
+DkJ-(f}+i/2) {-Zk kvi)] , 



where Ck and Dk are integration constants and where 
the complex number Zk is defined by 



- /I 



1^ 
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I arctan ^ 27/ fc^ ) 
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Eq. (154) should be compared to its non-CSL counter- 
part, Eq. (45). The only difference is the appearance 
of the Zk factor. This is consistent with the remark 
made above since this factor always multiplies the expres- 
sion krj and can, therefore, be viewed as a "renormaliza- 
tion" of the wavenumber k. In the non-CSL case where 
7 = 0, one obviously has zj, = 1 and Eq. (154) reduces 
to Eq. (45). It is interesting to remark that Zk for the 
parametric oscillator plays a role similar to that of a for 
the harmonic oscillator, see the definition (125). In fact, 
strictly following this last definition, one can introduce a 
mode-dependent ak parameter, namely ak = •\/2z7 — fc^ 
(using UJ = k for massless perturbations) and, then, Zk 
appears to be just a rescaled ak parameter: ak = ikzk- 
Finally, notice also that the sign ambiguity in the def- 
inition of Zk due to the presence of a square root has 
absolutely no impact on the results presented below. 

Let us now discuss the solution fk{'n) and what this 
implies for the behavior of the wavefunction. In pres- 
ence of the CSL term, the problem is characterized by 
three scales: the wavelength of the Fourier mode given 
by ^kiv) — '^{v)/k, the Hubble radius in (r/) = a^/a' and 
a new scale associated with the parameter 7 defined by 
£^ = a {rj) I ^ or, in terms of mass scale, = ^jaij]). 
Notice that is a new, time-dependent, physical scale 
that is built in the inflationary CSL theory, see Fig. 4. 
In terms of these three physical scales, the quantity Zkkr\ 
which appears in Eq. (154) can be written as 



Zkk^ = (1 + ,9) 



k"^ 
'^phys 



(156) 



where fcphys — kjavs, the physical wavenumber. At the 
beginning of inflation, the modes of cosmological inter- 
est today laid far inside the Hubble radius, which means 
Afc <C ^H, J-c. kri —00. Notice that these consider- 
ations are independent of the value of AI^. Indeed, if 
fcphys 3> M^, then Zk — 1 and the previous limit is not 
changed. On the contrary, if fcphys <C Mj, then the condi- 
tion \zk \ ^ 1 is even better satisfied. It is also interesting 
to remark that, in this last case, Zkkrj does not go to —00 
along the real axis but along a direction that is inclined 
in the complex plane. However, this does not change 
the asymptotic behavior of the Bessel functions in this 
regime. Upon using Eq. (154), one obtains 



lim fk{r]) 



Ck sin (^-Zk krj - ^/?^ 



+Dk cos {-Zk kri + 



(157) 



(154) This expression can also be re-expressed in term of "plane 
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wave" functions (writing = + ioi\.) 

lim fkiv) = 4^e"^l''l-'"'^"-*"/4 
Afc/£H->o V27r 

\/27r 

where the coefficients and can be expressed as 
linear combinations of Ck and Df., namely 

Ak = Cfce-^'^(^+i/2)/2 + i^^e"('3+i/2)/2 (159) 
Bk = Cfce"(^+i/2)/2 + 25^e-"(^+i/2)/2. (160) 

The solution (f58) is nothing but the Wentzel-Kramers- 
Brillouin (WKB) mode function exp(±i f ojdr) / ^/2uj. 
The reason for this result is that, in the sub-Hubble 
regime, the WKB approximation is still valid even in 
presence of the CSL term. As is well known, this ap- 
proximation is satisfied when the quantity \QluP'\ <^ 1, 
where Q is given by 

Since, in the limit under consideration, uP' tends toward 
a constant, namely uP' — k'^ — 2^7, and since Q is given in 
terms of derivatives of w, it is obvious that the criterion 
is satisfied. As already mentioned, the only effect of the 
CSL theory is to add the constant term — 217 to uP . Al- 
though this modifies the solution for the mode function, 
clearly, this cannot change the fact that WKB is valid at 
the beginning of inflation. 

Let us now comment on Eq. (158). When jryj goes to 
infinity, the second branch of the above solution is going 
to die away since > 0. As a consequence, only the 
first branch remains and, since $7^ is given in terms of a 
ratio, i.e. —if^/i^fk), the remaining constant Ap. disap- 
pears from the final expression. Therefore, fife becomes 
independent of the initial conditions and is given by 
rife ~ iak/2, which implies that Sftefifc ~ —a\/2 ~ —k/2. 
Returning to Eq. (f43), this means that the wavefunc- 
tion is not bounded at infinity and is not normalizable. 
The deep reason is that, in the CSL context, (or a^) 



is complex and this implies that the WKB solution ac- 
quires either a growing or a decaying exponential compo- 
nent which automatically kills one of the two branches. 
And, of course, (or ak) is complex because of the CSL 
term —2ij. 

Based on the previous discussion, it is clear that the 
only meaningful choice of initial conditions is to require 
that Ak — 0. From Eqs. (f 59) and (f60), we see that this 
implies 

Ck = -Dke''''-^+^P\ (162) 

This choice exactly coincides with the Bunch-Davies ini- 
tial conditions (46). From now on, we assume Eq. (f62) 
but we will come back soon to this discussion. Then, 
one can re-derive the behavior of in the sub-Hubble 
regime. One obtains 

1 

lim Slfc(77) = --Qffc, (163) 

Afe/^H^O 2 

which is fully consistent with Eqs. (f33) and (f34). In 
particular, one can check that, now, Sftefifc — )• k/2 and 
the wavefunction becomes normalizable (of course, it 
tends to the ground state wavefunction). Therefore, we 
have proven that, as expected, the cosmological pertur- 
bations behave, in the sub-Hubble regime, exactly as the 
CSL harmonic oscillator. 

Having studied the behavior of the stochastic wave- 
function in the sub-Hubble regime, we now turn to the 
super-Hubble case. In the framework of CSL, and con- 
trary to the sub-Hubble regime studied before, it is clear 
that this regime has no counter part in the case of the 
harmonic oscillator. It corresponds to the limit £h ^ Afc 
and, from Eq. (156), we see that this means |2:fcfc?7| — 0. 
Let us notice that one could also consider the case where 
fcphys ^ My such that A/^/fcphyg :s> 1 compensates the 
ratio iii/Xk in Eq. (156) resulting in a large Izfcfcryl, even 
in the super-Hubble regime. Below, we briefly comment 
on this case. Here, we assume that is such that this 
does not happen. Then, upon using the asymptotic be- 
havior of the Bessel functions for small values of their 
argument, one arrives at 



^(l + /3) ^{-kr,) j-kr,) j ( _ „,^\-^-V2 l\ TjP + 1/2) , .-2p-2, 

k~ 2krj 4(/3 + 3/2) 2{p + 3/2) k^^ Ck \ k^J ^ ^ 2 ) T{-p - 1/2)^ ^' ' ' 

(164) 

This equation should be compared to the corresponding non-CSL formula (52). If 7 = and if one takes the Bunch- 
Davies initial conditions, = —Cke~^^^-^~^^^^\ then the above equation exactly reduces to Eq. (52). Here, although 
we argued before that one should use the Bunch-Davies initial conditions (162), we temporarily keep the coefficients 
Ck and Dk arbitrary because, later on, we shall want to comment on their influence on the shape of the CSL power 
spectrum. Let us also notice that the last term of the above expression is in fact proportional to Zj^ (2/3-(-i)^ j£ write 
Zk in polar form, Zk = |2fc|e'^'= (of course, 9k should not be confused with the squeezing angle) where the modulus 
and the phase can be read off directly from Eq. (155), and parametrize the initial conditions as Ck = |Cfc|e'^'' and 
Dk = \Dk\e^^''~^'^^~^^'^P (so that the Bunch-Davies limit is simply 9d — dc = 0), then it is easy to determine the real 
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and imaginary parts of the function Qk. One finds 



Sm fife (77) 



T2(-^- 1/2) cos(7r/3) 



-2^-2 



X i sin + {2P + 1)6^-9^ + e,] cos(7r^)r2 (^/3 + (-fc,y)-2^"2 ^ 
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r 



These equations are the CSL counterparts of Eqs. (53) 
and (54). Of course, for 7 = and the Bunch-Davies ini- 
tial conditions, they exactly reduce to those equations. 
We see that the main effect of the CSL theory is to 
strongly modify 5Re since its leading term in the above 
expansion is a term which cancels if 7 = 0. We also see 
that we still have Jftefifc — >■ in the super-Hubble limit. 
In absence of the CSL term, we would obtain the same 
limit but not with the same power. Compared to ^eftk, 
^JmOfc is much less modified since the first correction 
show up only in the third term of the expansion. As 
a consequence, we still have Qm Cl^. — > 00 in the super- 
Hubble regime. 

We now use the above results to discuss the collapse of 
the wavefunction in more detail. Since we have assumed 
in Eq. (142) that the "collapse operator" is Vk, we ex- 
pect the non linear and stochastic terms in the modified 
Schrodinger equation to drive the initial Gaussian state 
to an eigenvector of v^, that is to say to the Dirac func- 
tion S (ufc — Vk). However, in practice, as we learned from 
the harmonic oscillator example in Sec. IV B, this is not 
what happens. In practice, we find that the wavefunction 
tends towards a Gaussian state with a constant spread 
in position and that the larger the value of 7, the smaller 
the amplitude of this spread, i.e. lim — > [h/ {4mj)]^^'^ 

when 7 — >■ CXI. Therefore, strictly speaking, the exact 
localization is obtained only in the 7 — >■ cxd limit. Of 
course, if the spread is very small, then for all practi- 
cal purposes, the collapse has been achieved. In fact, 
this is the essence of the amplification mechanism dis- 
cussed in Sec. IV A. The effective value 7^^ of 7 for a 
macroscopic object (or for its center of mass) is the fun- 
damental 7 times the number of particles in that object 
which results in a huge effective 7 and, therefore, a very 
efficient localization. As a consequence, a collapse can 
occur for macroscopic objects while it does not happen 
for microscopic particles even if their behavior is slightly 
disturbed. 

Let us now see how the previous discussion applies to 
inflation. The first difference is that the standard devia- 
tion, 1/(2-^3^6 rifc), does not go to a constant as for the 
harmonic oscillator but to infinity since Eq. (165) implies 
that rife cx 77 — )• 0. We remark that the divergence is less 
violent than when 7 7^ since, in that case, flk oc ry^ — >■ 0, 



according to Eq. (53). This is of course due to the in- 
fluence of the non linear and stochastic terms. However, 
this influence is not sufficient to prevent the divergence 
of the variance and, therefore, to ensure an efficient lo- 
calization. As a matter of fact, we see that, in the limit 
?7 — > 0, the main divergence in the Hamiltonian comes 
from the term ex w^ti^ while the CSL term goes like 7ti|. 
Hence, it is because the term c>c 77"^ diverges at the 
end of inflation that the Hamiltonian strongly dominates 
the dynamics of the system, preventing the CSL terms 
cx 7w| to carry out its job and to localize Vk (however, 
see Appendix A). This is certainly a problem for the in- 
flationary CSL theory. This issue can also be related to 
the fact that it is unclear how an amplification mecha- 
nism could be implemented in quantum field theory. As 
a consequence, the collapse mechanism is controlled by 
the parameter 7 and no effective 7 can be derived which 
would ensure a better localization. 

Finally, let us mention that one could wonder whether 
the localization can be achieved during the radiation 
dominated era that takes place after infiation. In this 
case, the scale factor behaves as 0(77) cx and, there- 
fore, ei = 2 and {ay/el)" / {ay/el) = 0. As consequence, 
the mode equation for fk is exactly that of an harmonic 
oscillator. This means that the variance now goes to a 
constant, see Sec. IV B, which seems to cure the problem 
discussed above. However, one can show that the corre- 
sponding value remains large for modes of astrophysical 
interest today. Therefore, this remains an unsatisfactory 
solution. 



D. The CSL Power Spectrum 

We now turn to one of the main goal of the present 
paper, namely the determination of the power spectrum 
predicted by the CSL theory. It was shown in Eqs. (29) 
and (37) that the power spectrum of the conserved quan- 
tity Cfe can be expressed as 



Vcik) = 



fc3 



1 



(167) 



Since we have determined the quantity 3f?e Qk in 
Eq. (165), the calculation of 7^^ becomes straightforward. 
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One obtains 



1 



Xn(fc)Ltand 



/3 + 3/2 



(168) 



where '^cl stand standard power spectrum given 

by Eq. (47) and the function /(/?) has been defined in 
Eq. (48). The function g-yik,^) 



cos(7r/3) 



cos [7r/3 + (2/3 + l)0k-Od + 9,] 

is seen to depend on the choice of the initial conditions. 
It has the property that, for 7 = and the Bunch-Davies 
initial conditions, gj=o{k,/3) = 1. In this case, and as 
expected, one can check that the modified power spec- 



trum (168) reduces to the standard inflationary power 
spectrum. We also notice that the power spectrum (168) 
is still a time-dependent quantity, contrary to the conven- 
tional case where the time dependence cancels out. For 
this reason, it is convenient to evaluate it at the end of 
inflation. In that case, the quantity —kr] can be rewritten 
as 



fcr, = --(l 
ko 



(170) 



where /cq is the comoving wavenumber of the Fourier 
mode, the wavelength of which equals the Hubble radius 
today, i.e. fco = ciqHq. The quantity AA^* denotes the 
number of e-folds spent by a mode of cosmological rele- 
vance today outside the Hubble radius during inflation; 
typically, one has AiV* ~ 50 — 60. As a consequence, the 
power spectrum (168) can be re-expressed as 



Vcik) = g^{k,p) 



J^g^{k,l3)f{l3) 



+ 



213+3 



(/3 + 3/2) 



.g(2/3+3)AiV./(l+/3) 



2/3+1' 



H -1 



Vcik) 



stand ' 



(171) 



Let us notice that, in Eq. (169), the quantities \zk\ of 
Eq. (155) and 9^ must now be written as 



1 



1/4 



-arctan 
2 



2^ f- 
kg \ k 



(172) 
(173) 



such that the amplitude of the CSL correction is con- 
trolled by the dimensionless ratio 7/^0- The for- 
mula (171) is one of the main results of this article and 
the corresponding power spectra for different values of 
the ratio 7/^9 are represented in Fig. 5. 

Let us now discuss in more detail the CSL power 
spectrum (171). Firstly, we notice that, in the short- 
wavelength regime fc/fco 00 j the power spectrum re- 
duces to V(;{k) ~ gj{k, l3)V(^\stand- Morcovcr, in this 
limit, we see that \zk\ 1 and 6*^ — 0. As a conse- 
quence, an almost scale invariant (namely, =2/3 + 4: 
with /3 < —2) power spectrum is recovered if one as- 
sumes the Bunch-Davies initial conditions, |Cfc| = |Dfc| 
and 9d — 9c = since, in that case, g^{k,(3) = 1. This 
almost scale invariant branch of the power spectrum is 
clearly seen in Fig. 5. Secondly, there is clearly another 
regime which corresponds to the case where the second 
term in the square brackets in Eq. (171) starts playing a 
role. If we neglect factors of order one, this happens at 



kj, where k^ solves 



^g,(fc„/3)e(^^+3)^^-/(i+«''^-' 



2,3+1 



(174) 



The value of g^ is mainly controlled by the value of \zk\ 
which is always close to unity provided that k <^ k^ with 



ko 



V2 



1/2 



(175) 



Then, let us assume that g^ — 1 when the condition (174) 
is met. In this case, the scale k^ can be expressed as 



"■7 

ko 



-1/(1+2/3) 



exp 



2/3-1-3 
-l)(2/3- 



1) 



-AN, 



(176) 

Choosing the fiducial value /3 ~ — 2 leads to k^/ko ^ 
(7/fco)^/'^ exp (AA^*/3). One can check that, indeed, 
kj ^ kz and, therefore, assuming g-^ ~ 1 was, in retro- 
spect, valid. As a consequence, in the range k <^ k^, the 
spectrum approximately behaves as oc = 
k^, that is to say with a spectral index of rig ~ 4. This 
second branch is also clearly visible in Fig. 5. In addi- 
tion, the dependence in g^ is canceled out which means 
that this prediction is actually independent of the choice 
of the initial conditions, a remarkable property indeed 
(this also means that, even if fc ^ fc^, the spectral index 
remains the same). Moreover, we see that this spectral 
index is also independent of /3 which is also remarkable. 
In this sense, the CSL branch of the power spectrum can 
be said to be "universal" (unfortunately not scale invari- 
ant!). 

We are now in a position where we can discuss the cos- 
mological constraints on the parameter 7. From the high 
accuracy measurements of the CMB anisotropics [20-22], 
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FIG. 5: Ratio of the power spectrum given by Eq. (171) to the standard power spectrum given by Eq. (47) for different values 
of the parameter 7/fco (and for /3 — —2.01, a value leading to a standard power spectrum close to scale invariance). The number 
of e-folds between Hubble radius crossing and the end of inflation (for the modes of cosmological interest today) has been taken 
to AA^t = 60 and the initial conditions have been chosen to be the adiabatic vacuum. 



we know that the power spectrum is almost scale invari- 
ant, ~ 1, and that a spectral index rig = 4 is com- 
pletely excluded. This means that the CSL branch must 
correspond to scales much larger than the present Hubble 
radius, in other words k^/k^ <C 1. This condition means 
that, for /? ~ — 2, one has 



constraint (177) can be rewritten as 



1_ 

/4q 



< e 



'AN, 



10~ 



(177) 



To our knowledge, this is the first time that a constraint 
on the parameter 7 is obtained from cosmological con- 
siderations (see, however, Ref. [TG]). We see that the 
constraint is expressed as a limit not on 7 itself but on 
the combination 7/fcg where we remind that fco is the co- 
moving wavenumber of the Hubble radius today. Looking 
at Eqs. (142) and (152), this was expected since the CSL 
modification amounts to a redefinition of the comoving 
wavenumber fc^ — >■ fc^ — 217. This means that, in order to 
characterize the amplitude of the modification, one has 
to compare the comoving wavenumber to 7, hence the 
ratio 7/fcQ. The appearance of the comoving wavenum- 
ber in the observational constraint reflects the fact that 
the theory contains a built-in "time-dependent physical 
preferred scale" i^yif]). In terms of physical scales, the 



< 10 



-13 



(178) 



today 



Clearly, the constraint is very strong and means that the 
scale £j is very large in comparison to the Hubble ra- 
dius today. This is another illustration of the fact that 
squeezed states are fragile and easily perturbed. For the 
CSL theory itself, this probably means that, in order 
to be compatible with cosmological inflation, an impor- 
tant fine-tuning is required. Of course, this conclusion 
should be toned down given the uncertainty that exists 
on a CSL formulation of quantum field theory as dis- 
cussed in Sec. VA. One might argue for instance that 
the above result could be due to the fact that our mod- 
ified Schrodinger equation is not necessarily the appro- 
priate one in the context of quantum field theory. It 
would also be interesting to compare the cosmological 
constraint with the other constraints on 7 derived in 
the literature. But, as explained before, because we as- 
sumed Vk to be the preferred basis for the collapse, our 
parameter 7 is actually different from the parameter 7 
considered elsewhere, in particular it has a different di- 
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mension. This complicates tremendously any comparison 
with other systems. 

Finally, before closing this section, let us discuss the 
following question. In this article, we have defined the 
power spectrum in the CSL theory by means of the for- 
mula E — E((?)fc)^). However, there is an issue 
regarding this definition. Indeed, it is clear that it does 
not go to zero when the parameter 7 vanishes. Ac- 
tually, it tends towards the standard result when the 
Schrodinger equation is recovered. However, it was ar- 
gued in Ref. [112] that the power spectrum should go 
to zero in the limit where 7 — >■ and, therefore, can- 
not be given by the definition used above. The reason 
advocated by Ref. [112] is that, without a collapse, the 
theory remains homogeneous and isotropic and, as con- 
sequence, there is simply no perturbations at all. This 
has led Ref. [112] to define the CSL power spectrum by 

E (^{vk)^^ — {{vk)), a quantity which indeed vanishes 
when 7 — and differs from the previous one. In this 
last paragraph, we explore the difference between these 
two alternative definitions. 

At any time, the wavefunction can always be expanded 

as 

* iv, Vk) = y" * iv, «fc) S (vk - Vk) dvk , (179) 

where the superscripts "R, I" have been ignored for con- 
venience. If a dynamical collapse of the wavefunction 
takes place then ^ is projected (collapsed) on an eigen- 
statc of the operator Vk, namely 

* ^ *eol = S{vk- Vk) , (180) 

where Vk depends on the specific realization under con- 
sideration, then one obviously has 

(^'col «'col) = (*col|wfc|lAcol)^ 



Therefore, for each realization, one has (wl) — (w/c)^, 
once the wavefunction has collapsed. Since this is true 



for all realizations, it remains the case after taking the 
stochastic average. Therefore, after the collapse, one can 
write 

E{{vl))=E[{vk)') , (182) 

and this remains true for any Hermitian operator. Note 
that this argument strongly depends on the fact that the 
wavefunction has actually collapsed to an eigenstate of 
the operator Vk- For instance, in the case of an harmonic 
oscillator studied in Sec. IV B, it was shown that the 
asymptotic state is not exactly a Dirac wavefunction, but 
a Gaussian state the spread of which does not vanish 
for finite values of 7. In that situation, the two above 
expressions are not identical. 

On the other hand, the second terms in both definitions 
of the power spectrum differ 

E(^{vk)^^ ^E^ivk)) , (183) 

so the two spectra do not coincide even after the col- 
lapse. The difference ultimately boils down to the fact 
that it is built out of a standard deviation which is not 
an Hermitian operator. This is a generic question for the 
predictions of any theory mixing different kinds of aver- 
ages (in the case at hand, quantum and stochastic) when- 
ever non-linear combinations of Hermitian operators are 
involved. 



VI. THE COLLAPSE OF COSMOLOGICAL 
PERTURBATIONS 

In this section, we investigate the collapse mechanism 
and its dynamics in more detail. In particular, we calcu- 
late the collapse time and compare it with the cosmolog- 
ical characteristic times. For this purpose, we now con- 
sider the following double Gaussian quantum state [73] 



^k{v,vk) = \Ni'Uv)\exp{~^eni'Uv) 



Vk - vl^^ (ry) 



.(1) 



iv) + 'iXk\v)vk 



A^f ^ (77) I exp{-3?er!if^ (77) [vk - «f (v)] + i<ji^'^ (v) + iXk\v)vk - t^mn^^\fj) (w^)'}, (184) 



where, as before, [A^^^'^'l, 3fief^^^'^\ ^i^'^^ '^k'^^ ^ Xk'"^^ ^'^'^ 3m 17^^'^^ are real, possibly stochastic, numbers. The 
superscripts "R, I" have not been written for convenience but it should be remembered that they are of course present. 
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Inserting the above state into the modified Schrodinger equation leads to the following set of formulas 



(1:2) 



(1,2) 



5m 



(1,2) 



.(1,2) 



Jl,2) 



(1,2) 



(1,2) 



d?/ 



(1,2) 



7 + 4 



3fierj 



(1,2) 



Sm 17 



(1,2) 



(1,2) 



Sm 



(1,2) 



(1,2) 



V7 dW^, 



5m r2 



(1,2) 



-5Rel7 



(1,2) 



= -4 5Rer2 



)(1,2) 



SRer? 



,(1,2) 



(1,2) 



2 r 



,(1,2) 



1 r 

2 



(1,2) 



(1,2) 



(1.2) 



2x 



(1,2) 0^f^(l,2) 



(1,2) 



(185) 

(186) 
(187) 
(188) 

(189) 
(190) 



These equations should be compared to Eqs. (144), (145), 
(146), (147), (148) and (149). They are obviously very 
similar except the two last terms of Eq. (185) and the 
last term of Eq. (188) which are new. In the case of a 
single Gaussian, one has (vk) = Vk and these terms disap- 
pear. In the present case, the expression of {vk) is a very 
complicated function of all the parameters describing the 
wavefunction. Let us also notice that, since the evolution 

of CTfc^'^^ and Xfc^'^^ depends on v\^''^\ these quantities 
also feel the coupling between the two Gaussian compo- 
nents. However, one can see that the equations of mo- 



tion for Jftefi^^'^'* and 3m $7^^'^'' decouple from the other 
equations of motion and form an independent and closed 
subsystem. This means that the evolution of these two 
functions is identical to that of their counterpart in the 
simple Gaussian case and, moreover, that, if the initial 
conditions are chosen to be the same, ^^j^^ — fi^^' at any 
subsequent time. From now on, for this reason, the su- 
perscripts "(1)" and/or "(2)" on these quantities will be 
dropped. 



(1,2) 



It should be clear that the above system of differential 
equations is rather complicated to study. However, as 
we shall see, the most relevant properties of the evolu- 
tion of the double Gaussian quantum state can be ana- 
lyzed in a rigorous way. In particular, it is interesting to 



In 



(2) 



N, 



(1) 



introduce the function Tk{i]) 

Ref. [73] . This quantity characterizes the relative impor- 
tance of one Gaussian component to the other and, there- 
fore, provides a criterion to decide whether the collapse 
has taken place. The superposition of the two Gaussian 
quantum states reduces to one of them when | Vk \ goes 
to infinity. In practice, the collapse will be said to have 
occurred when | \> h with, say, 6 ~ 10 [73]. Then, by 
subtracting the two equations (185), one arrives at the 



following evolution equation for F^ 



dFfc 
d?7 



-(2) -(1) 
^k -^k 



7 



-(2) -(1) 
^k -^k 



2{vk) 



.(191) 



This equation remains complicated because of the pres- 
ence of the term (ufc). However, the calculation 
can be simplified if one assumes that the two Gaus- 
sian components of the wave function do not overlap, 
i.e. have separate supports. Technically, this means that 



SRerjfe 
formula 



-(2) -(1) 



S> 1, leading to the following simple 



{Vk 



(1) 



,(1) 



(2) 



,(2) 



N, 



(1) 



(2) 



(192) 



Inserting this formula into Eq. (191) and defining Xk by 



Xk = V 



- „-,(2) 



vl^\ one obtains the following expression 



^ = VlXk^+lXltanh{Tk) . (193) 

This stochastic differential equation can be further sim- 
plified. Indeed, using the new timelike variable [73] 



Sk = l I Xl (u) du , 
'mm 



Eq. (193) can be rewritten as 
dFfc dWs , 



dsk dsk 



tanh (Ffe) 



where 



Ws=Vl I Xk dWr, 
'0 



(194) 



(195) 



(196) 



is another Wiener process with respect to the time vari- 
able Sfc. 
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A. Collapse Time: Definition 

Let us now study the stochastic differential equation 
driving the evolution of Tk in more detail. In particular, 
we would like to know how much time it takes for the 
wavefunction to collapse or, in technical terms, we would 
like to determine the value of Sk such that Ir^l > b. The 
quantity Tp, being stochastic, two complications arise. 
Firstly, once it has reached a value larger than b, there 
is no guarantee that it will stay in this region. The ran- 
dom behavior of Ffc could temporally brings it back to 
the region \Tk\ < b. However, since the average trend is 
clearly to have a collapse, this would happen for a lim- 
ited amount of time only before returns in the regime 
where Ir^l > b. For this reason, we will consider that the 
wavefunction has collapsed when F^ has crossed the value 
±6 for the first time. Technically, this means that we are 
led to define the "collapse time", Sk, as Sk = inf (sfc) 
such that |Ffc {sk)\ > b, see also Ref. [73]. A second issue 
is that, clearly, the value of Sk will differ from one re- 
alization to the other or, in other words, that Sk is still 
a random variable. Therefore, we will rather define the 
collapse time as the ensemble average value of Sk but 
we will also be interested in calculating its higher order 
momenta. 

We now seek an explicit expression for the quantity 
Sk- It can be obtained in the following manner. Let 
us consider a function c(Ffc) that we do not characterize 
in more detail for the moment (but see below). It can 
always be Taylor expanded in dF^. At second order, the 
result reads 

c(Ffc + dFfc) = c(Ffc)+c'(Ffc)dFfc 

+ lc"iTk)drl + 0{drl) , (197) 

where dFfc is given by Eq. (195). At first order in dsk, 
this leads to 

dc[Ffe(sfc)] = c'[Tkisk)]dWs 

+ c' [Tk{sk)]ta.nh[Tk{sk)]dsk 

+ ^c"[Ffc(sfe)]dsfc. (198) 

Then, integrating the above expression between Sk = 
where F^ {sk = 0) — 69 and Sk — Sk where 
Ffc (sfc = Sk) = ±&, one gets the following (Ito) formula 

c{±b)-c{bo) = / c'[rk{sk)]dW, 
Jo 

(199) 

At this stage, we now specify the function c. We require 
it to be the solution of the differential ordinary equation 

^c"(a;) -)-tanh(a;)c'(a;) = -1, (200) 



with boundary conditions c {—b) = c {+b) = 0. It is easy 
to show that c(x) = 6tanh(6) — a;tanh(x). This means 
that the first term on the left hand side of Eq. (199) 
vanishes and that the integrand of the second term on 
the right hand side is just —1. Therefore, Eq. (199) can 
be rewritten as 

Sk = c (bo) + f " c' [Tk (sk)] dW, , (201) 

JQ 

and this gives an (implicit) expression for the quantity 
Sk- Finally, by averaging over all realizations, one ob- 
tains [73] 

E (Sk) = c (bo) = 5tanh(6) - 60 tanh(6o) • (202) 

The fact that the stochastic average of the integral in 
Eq. (201) vanishes comes from the fact that c' [F^ (sfc)] 
depends only on stochastic events occurring at s'^ < Sk- 
As a consequence, it can be expressed as an integration 
over ds^ and dWs' where s'^ < Sk- Since E (dW^^'dM^s) = 
S (sj, — Sk) ds|, at first order in dsk, the stochastic av- 
erage of the integral term in Eq. (201) vanishes. Actu- 
ally, things are slightly more complicated since the upper 
bound of this integral, S'fc, is a stochastic quantity itself. 
Therefore, the averaging process should also be carried 
out on this upper bound, and a generalized demonstra- 
tion which includes this case can be found in Ref. ([113]) 
(theorem 1 on p. 28). 

In order to characterize better the properties of this 
collapse time, it is also important to determine its vari- 
ance. Interestingly enough, the same technique described 
above can be used in order to calculate iteratively higher 
orders of Sk- Upon using Eq. (201) one has 

rSh 

E (si) - (bo) + / c'2 [Ffc (sk)] dsk - (203) 
Jo 

We see that we now need to evaluate the integral in the 
above expression. For this purpose, we consider a new 
function e(Ffc). As was done before, it can be Taylor 
expanded and this leads exactly to Eq. (199) (with, of 
course, c replaced by e). Compared with the proof that 
allowed us to obtain E(S'fc), at this point, the strategy 
changes. We now require the function e (x) to be the 
solution of the following ordinary differential equation 
[compare with Eq. (200)] 

^e" {x) + tanh (x) e' (x) = -e'^ (x) , (204) 

with boundary conditions e (—6) = e{b) = 0. As before, 
one can use this differential equation into the Ito formula 
to simplify the second integral in Eq. (203) [more pre- 
cisely, the integrand is replaced by — e'^ (x)]- Taking the 
stochastic average of the resulting equation, one gets 

e (60) = / ' e'2 [Ffe (sfc)] dsk - (205) 
Jo 
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As a consequence, we deduce that 

E (Si) = (60) + e (bo) 



(206) 



The only thing which remains to be done is to solve 
Eq. (204). In fact, it turns out to be more convenient to 
solve the slightly simpler differential equation satisfied by 
d (x) = {x) + e (x), namely e"(a;)/2-|-tanh (x) e'l (x) = 
—2c{x), with boundary conditions ei (—6) = ei (b) = 0. 
It is straightforward to show that 61(2:) ~ — 6^ + 
[1 + 26 tanh (b)] [b tanh [b) - x tanh (x)]. Then, the sec- 
ond momentum of S can be simply expressed as E {S^) — 
ei (60) which, therefore, gives an explicit expression for 
the variance of the collapse time. Since b is supposed to 
be a large number b ^ 1 and if we assume that the two 
Gaussians have comparable initial weights which implies 
that bo ^ 0, then one obtains, at leading order in 6, 



E{Sl)~E^Sk) ^ Vb. 



(207) 
(208) 



These two equations tell us that the relative standard 
deviation scales as l/Vb and, therefore, that the distri- 
bution of Sk becomes more peaked as b increases. For 
this reason, in the following, we will simply estimate the 
collapse time by means of the sloppy requirement that 
Sk = b. Finally, let us mention that one could also apply 
the technique used in this section in order to determine 
the higher order correlation functions of the process Sk. 



B. Collapse Time in the sub-Hubble Regime 

In the last section, we have explained how to determine 
the collapse time in terms of the variable s^. In order 
to translate this result in terms of a more physical time 
(conformal time or, better, number of e-folds), we need 



to use Eq. (194) which, in turn, requires the knowledge of 
the function Xk- This one cannot be determined in full 
generality but it is easy to characterize it in the sub- and 
super-Hubble regimes. In this section, we investigate the 
sub-Hubble regime. 

Let us define Kk = xi^"* ~ Xk'^ ■ This quantity mea- 
sures the shift in momentum between the two Gaussian 
components of the wavefunction (184) (we recall that Xk 
measures the shift in position). Then, taking the differ- 
ence between the versions "(1)" and (2)" of Eq. (188) on 
the one hand, and versions "(1)" and (2)" of Eq. (190) 
on the other hand, we arrive at a closed system which 
can be written in a matrix form, namely 



d 
drj 



Xk 
Kk 



-23ni Qk 



7 



1 



^e^k 

Ai^erikf 25mr2fc 




(209) 

At this stage, there is no approximation and the above 
equation is general. In the sub-Hubble regime, one can 
use Eq. (163) to simplify the expressions of dieftk and 
Smfifc. Moreover, we are mainly interested in com- 
puting the collapse time for the modes that correspond 
to the (almost) scale invariant part of the power spec- 
trum since it is clearly less interesting to compute this 
quantity in a regime that is already excluded by the 
data. As was discussed before, this amounts to consid- 
ering that ^ 1. Under those conditions, one has 
^eflk -> k/2 and 3m fife -7/(2fc) and Eq. (209) can 
be re-expressed as 



±f Xk 
drj \ Kk 



-7/fc 1 



Xk 
Kk 



(210) 



This system of differential equations can be integrated 
and the solution reads 



Kk iv) = e-^(''-''"")/'=<^ Kk,ini cos [fc (r; - yyini)] - kXk^i sin [fc (r; - yyini)] 



(211) 



Xk iv) = e-^("-""")/'=<{ XkMi cos [fc - mm)] + sin [fc {rj - mni)] 



(212) 



where Kk^ini and Xk^ini are two integration constants conveniently chosen to be the values of Kk and Xk at initial 
time 77 = rjini- For simplicity, we now consider a situation such that Kk.mi = 0. Upon using Eq. (194), one finds that 



Sk 



:Xt 



-27(r;-jji„i)/fe 



1 



7 ^2 
fc3 



1 



'1-1-472/^4 



-27(r,-r,i„,)/fc {cos [2fc (t; - r?i„i)] - siu [2fc (77 - 7?i„i)] - 1} . 

(213) 



If we expand the above result in 7/fc^ for the reason approximated expression for the mapping between the 
discussed before then, at leading order, one obtains an 
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variables ry and Sk 



Sk 



-27('i-')ini)/fc 



(214) 



This expression means that Sk runs from to kX^ ini/4 
when 77 runs from yyini to infinity. Therefore, the time 
Sk evolves in a finite range. However, in order to be 
consistent, one must have 77 < 77* = — 1/fc since the 
equations that have been used in order to derive Sk are 
vafid only in the sub-Hubble regime. As a consequence, 
we have in fact Sk £ [0, s*] where s* = kX^ ^^J4:{1 — 
exp[(27/A;^)(l + kijini)]}- Since we have |fc77ini| ^ 1, 
one can thus write s* ~ kX'^ ini/4[l ~ exp(2777ini/fc)]. If 
s < s*, then Eq. (214) can be inverted in order to eval- 
uate the (total) number of e-folds in terms of the time 
variable Sk- One finds 



Nk = (l + /3)ln 



1 - 



1 



27 fc77ii; 



In 1 



4sfc 



kXi 



(215) 



and one checks that if Sfc = then Nk = 0, if Sk ~ s* 
then — > 00, and that the condition s < s* is sufficient 
to guarantee that the above expression is well defined. 

Let us now discuss the above results in more detail. 
Firstly, we notice in Eqs. (211) and (212) that the func- 
tions Kkiji) and Xk(j]) tend to zero when 77 — 77ini 3> 1. 
When this happens, the two Gaussians have the same 
mean in position and momentum; in other words the two 
Gaussians have merged. This "merging phenomenon" 
seems to be a generic feature and can also be observed 
for the free particle [73] and/or the harmonic oscillator 
in Minkowski spacetime. Therefore, it does not come as 
a surprise that it also shows up in the sub-Hubble regime 
where the Fourier mode under consideration does not feel 
spacetime curvature. This also means that it is not a pe- 
culiar property of inflation. 

The free particle situation can be studied [73] by 
returning to Eqs. (117), (118), (119), (120), (121) 
and (122). It is sufficient to consider that w = in 
those equations to obtain this case. This means that the 
mode equation (124) now reads — a^/fc — 0, where 
the quantity a, defined in Eq. (125) for the harmonic os- 
cillator, now reads a = ^J^iriKJm = y/jH/m{l + i) and 
is obtained from Eq. (125) by taking to = 0. As a conse- 
quence, the solution for n{t) has exactly the same form 
as in Eq. (126) but now with the new a given above. This 
implies that 3?eSl — > ^J^jmjhj'l and 9m O — >■ 
when t — > CX3. These formulas should be compared 
to Eqs. (133) and (134). Then, considering the equa- 
tions of motion for a double Gaussian state, and defining 
X = X2 — xi and K = X2 — Xi^ upon using Eq. (209), 
one obtains the following set of equations 



^( X 
At \ K 



-7 



m h/m 



(216) 



This equation should be compared to Eq. (210). In 
particular, one notices that, here, the free particle case 
is not simply obtained from this equation by consid- 
ering k — uj — If we assume that if(0) = 
0, then the solution for X(t) is given by X(t) = 
A(0) ex\){—tyjh'y I'm) cos{ty^h'y/m). We see that this so- 
lution resembles the solutions (212) and (211) obtained 
before. Therefore, the merging is indeed already present 
for a free particle in flat spacetime and is not a specific 
feature of inflation. The exponential factor is mainly re- 
sponsible for the merging and this means that the "merg- 
ing time" of the free particle is given by 



merge \l ■ 



(217) 



This expression is consistent with the merging time de- 
rived in Ref. [73]. 

In order to discuss our inflationary result, one should 
consider the merging time of the harmonic oscillator in- 
stead of that of the free particle since this is the appropri- 
ate limit in the sub-Hubble regime. Following the same 
logic as before, it is easy to show that, for the harmonic 
oscillator, Eq. (216) is replaced by 



d 
dt 



—^h/{muj) h/m 
—muP'/h —'-fh/{muj) 



(218) 

We see that it is indeed similar to Eq. (210) if we take 
u! — k (and m = h ^ 1). The solution for X{t) can be ex- 
pressed as X{t) = X{0) exp[—h'^/ {muj)t] cos(wi), assum- 
ing as before K{0) = 0. This solution is perfectly consis- 
tent with (211) and (212). Compared to the free particle 
case, one notices that the coefficient in the exponential 
is now different from the frequency of the trigonometric 
function. But the most important result that one can de- 
duce from the above considerations is that the merging 
phenomenon is also present for the harmonic oscillator 
and that the corresponding merging time is given by 



T 



ho 

merge 



V merge/ 



(219) 



Let us remark that the last expression could have been 
guessed on dimensional grounds. 

In the case of inflation, the conformal merging time is 
given by [see Eqs. (211) and (212)] 



k{^_ 



merge 



- ??ini) = 



fc2 

7 



(220) 



However, there is a new twist in the discussion. It is not 
obvious that the above equation admits a solution be- 
cause, in some sense, we have a limited amount of time 
from 77i,ii to 77* , the time of Hubble horizon crossing (de- 
fined by |fc77*| — 1). For times such that |A;77| < 1, we 
are no longer in the sub-Hubble regime and the above 
equation can no longer be used. But, given a value of 
fc^/7, and an initial time 77ini, it is not obvious that there 
exists a time 77meige such that Eq. (220) is satisfied. In 
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fact, there exists a solution only if |fc?7ini| > 1 + 
This condition means that, for a given one can al- 

ways give more time to the system to satisfy Eq. (220) 
by starting its evolution earlier (which is equivalent to 
increasing |77im|)- It is easy to show that the previous 
inequality is in fact a condition on the total number of 
e-folds during inflation (/? < —2), namely 



scales, the mode under consideration must behave as in 
flat spacetime. Indeed, for a free particle or the harmonic 
oscillator in Minkowski spacetime, the condition for the 
collapse to occur can be written as s = 7 / X^(r)dr ~ 
7X(0)2r*i'''° = b, where we have used X{t) ~ X{0) since 
we have shown that the merging takes place on a much 
longer time scale. This implies that 



Nt > AN^ + In 1 



7 



(221) 



where AA^* ~ 50 for the modes of cosmological interest 
today. If this condition is met, then the merging occurs 



after 7V™"«^'= with 



= - In 1 



(222) 



Moreover, the term k'^ / {'ykri^^i) is of the order ^ 
fc^e~^'^"'"^°/7 and it seems reasonable to assume that it 
is small. Indeed, typically, the total number of e-folds 
during inflation is very large and, even if fc^/7 3> 1, the 
factor e~^'^ will entirely compensate its influence (to be 
more concrete, we know that fc^/7 > 10^® but Nt can 
easily be larger than, say, 1000 and can even be as large 
as 10^). Then, the merging time during inflation can be 
approximated by 



jymorgo 



fc2 



< 1. 



We see that this expression scales as oc A;/7 in full agree- 
ment with the previous considerations on the harmonic 
oscillator, see Eq. (219). 

Let us now study the collapse time. First of all, the col- 
lapse can occur in the sub-Hubble regime only if 6 < s*. 
If we use the expression of and assume, as before, that 
k'^ / {jkrjini) <C 1, then ~ kX^^^jA and the condition 
for having the collapse in the sub-Hubble regime can be 
simply rewritten as 



Tl 



fp,ho 



col 



7^(0)2 



(226) 



and one verifies that it is similar to Eq. (225). Therefore, 
if the collapse occurs on sub-Hubble scales, its properties 
are, as expected, similar to what happens in flat space- 
time. Finally, if one starts from an initial state made of 
several well-separated Gaussian wavefunctions, the pre- 
vious calculation suggests that it will almost instanta- 
neously turn into a single Gaussian state. As a matter 
of fact, it is a general property [64, 114] of the CSL dy- 
namics that it asymptotically leads to Gaussian states. 
A posteriori, this remark reinforces the assumption of us- 
ing a Gaussian state for the calculation of the spectrum 
in Sec. VD. 

When the condition (224) is not satisfied, there will be 
no collapse on sub-Hubble scales. However, we can still 
hope it will happen on super-Hubble scales. In fact, the 
claim that the collapse has occurred depends on the value 
chosen for b. Before, we used 6 ~ 10 and for this value, 
(223) given that our working assumption is kX^ j^^j ^ 1, the 
condition (224) is probably always satisfied. Therefore, it 
is only if we are more demanding about the criterion that 
defines the collapse that this condition can be violated. It 
is clear that a more stringent criterion takes more time to 
be satisfied and, in this case, the time "at our disposal" 
in the sub-Hubble regime may not be sufficient. In this 
situation, we have to consider the super-Hubble regime. 
In the next section, we turn to this case and show that 
the collapse is less efficient on large scales than it is on 
small scales. 



6< 



kX'^ 



(224) 



If this condition is satisfied, then the "e-fold collapse 
number" of the mode under consideration is obtained 
by putting = 5 in the above expression (215). Upon 
using the same assumptions as before, we obtain that 



2b 



< 1. 



(225) 



At this point, several remarks are in order. Firstly, we 
notice that N^°^ /N^'''^'' = 46/(fcX| ;„;) < 1. This means 
that the collapse occurs on a much smaller time scale 
than the merging. This property was also noticed in the 
case of a free particle in Ref. [73]. This means that the 
merging cannot be viewed as a substitute for the collapse. 
Secondly, we notice that Nj^°^ is actually independent of 
k. We interpret this fact as meaning that, on sub-Hubble 



C. Collapse Time in the super-Hubble Regime 

In this section, we repeat the previous discussion but 
now in the super-Hubble regime. Therefore, we restart 
from the equations (209) but now use the super-Hubble 
limit (165) and (166) for ^eflk and SJmfifc. For the 
modes of cosmological interest today in the (almost) scale 
invariant branch of the CSL power spectrum, one has 
7/fc^ ^ 1 and the solution for Xk{ri) can be simply 
written as 



(227) 



where X^* is the value of Xk{ri) when the mode under 
consideration k crosses the Hubble radius. One can see 
that Afc (77) increases with time contrary to what happens 
in the sub-Hubble regime. From this expression, it is easy 
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to derive the relation between Sk and the conformal time. 
One obtains 



Sk = 



7 



A:2 2/3 + 3 



{-krif^+^ - 1 



(228) 



The last formula is valid only on super-Hubble time, 
that is to say for 77 > 77* = — 1/fc. At 77 = 77*, Sfc = and 
then Sfc — >■ 00 as 77 ^ 0. From this expression, it is also 
possible to relate the time variable Sk and the number of 
e-folds. One arrives at 



Nk^N, 



14 
2/3- 



-Sk 



(229) 



This expression is always well defined because 2/3 + 3 < 0. 
One verifies that Sfc = corresponds to Nk = iV* . 

Let us now derive the time of collapse. As usual, it is 
obtained by Sk = b. As a consequence, it is simply given 

by 



Nk^ 



N, 



1 + 13 
2/3 + 3" 



In 1 



fc2 2/3 + 3. 



(230) 



As a first check of this equation, we notice that, when 
7 — >■ cx), iVjg°' ~ A^*. Of course, this result is ex- 
pected since a large value of 7 means that the col- 
lapse mechanism is very efficient and, therefore, that 
the wavefunction almost instantaneously collapses. On 
the other hand, the formula (230) can be further simpli- 
fied. Indeed, if the collapse has not taken place on sub- 
Hubble scales, it is also the case for the merging since 
N^oi/j^mergc ^ ^ ^ couscqueuce, Xkiv) has not 

evolved much and one can replace Xk* by Xk.ini- More- 
over, for the same reason, one must have b > kX^ ini/^' 
see also Eq. (224). In addition, we know that k'^ /'j ^ 1. 
Therefore, the first term in the argument of the logarithm 
in Eq. (230) can be neglected. For /3 ~ —2, this equation 
can be rewritten as 



D. The Born Rule Derived 

Finally, we conclude with a section where we calculate 
the probabilities of collapsing to each of the two branches 
of the wavefunction. We show that these probabilities are 
given by the Born rule, which is of course expected since 
the CSL theory is precisely designed to reproduce this 
result, as already discussed in Sec. IV (see also Ref. [T-i]). 

Let us denote by pi the probability that the system 
collapses on the first Gaussian branch of the wavefunc- 
tion. This is also the probability that, from given initial 
conditions, the stochastic quantity Ffc reaches first the 
region F^ < — 6 (i.e. before the region F^ > 6) and that, 
therefore, one has Tk{Sk) = —b. Clearly, the probability 
P2 that the wavefunction collapses on the second branch 
is the probability that Ffc(5fc) = b. Now, let us introduce 
a function i}'{x) which is defined by 



g{x)-g{b) 

9{-b)~g{b) 



(232) 



where g{x) will be specified soon. By construction, one 
has 1/' {—b) = 1 and -0 (5) — 0. Since, by definition, 
Ffc(S'fc) can only take two values (namely ±6), one has 

E {7^ [Ffc {Sk)]} = Pi^ i-b) + P2^J (b) = pi, (233) 

and this gives us a method to calculate pi. To do so, we 
follow what was explained in Sec. VI A, see in particular 
Eq. (199), and we write the corresponding Ito formula 



i^[Tk{Sk)]-^{bo) 



^' [Tk (sfc)] dWs 



^' [Tk (sk)] tanh [Tk (sk)] + [Tk (sfc)]| dsfc. 

(234) 

Then, let us choose the function g(x) such that it obeys 
the equation 



1 



.g" {x) + tanh (x) g' (x) = , 



(235) 
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Of course the result will depend on what we require for 
b and what we assume for Xk.ini- However, it seems rea- 
sonable to assume that the second logarithm will not lead 
to a dominant contribution. If this is the case, then 
our result simply says that the wavefunction collapses 
just ln(fc2/7) e-folds after the Hubble radius crossing. 
Given the constraint obtained from the measurement of 
the power spectrum in Eq. (177), one already knows that 
N^°^ — > 28. Smaller values of 'y/k'^ would of course 
lead to a larger number of e-folds. We conclude this sec- 
tion by noticing that the constraint (177) is compatible 
with a collapse occurring during inflation. Only for val- 
ues of 7 such that 7/fc2 < 10-50 (and b > kXl .j4) 
would the collapse happen after inflation. 



or, equivalently, g{x) = tanh(a;). Since Eq. (232) im- 
plies that '4>{x) and g{x) are linearly related, '4>{x) also 
obeys the above differential equation. As a consequence, 
the second integral in Eq. (234) vanishes. Taking the 
stochastic average, one obtains 



E{7^[Ffc(5fc)]}=Pi =V(6o) , 



(236) 



which is explicitly known since g{x) has been determined. 

The probability p2 can be deduced along the same 
lines, by introducing a new function -0 such that, this 
time, ip{—b) — and ip{b) = 1. Another method, much 
simpler, is just to use the condition pi + P2 = 1- The 
flnal result reads 



Pi 



P2 



tanh (bo) — tanh (6) 
tanh (—6) — tanh (b) ' 
tanh (60) — tanh (— &) 
tanh (5) — tanh (—6) 



(237) 
(238) 
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From the definition of Ffc, these two formula can be 
rewritten as [73] 



Pi = 

P2 = 



\Ni (77i„i) 



\N, (77i„i) P + \N2 iVini) P ' 
|iVl (^ini) P + 1^2 (^ini) P ' 



(239) 
(240) 



which are exactly the Born rules of conventional quantum 
mechanics. 



VII. CONCLUSION 

Let us now summarize our main findings. In this pa- 
per, we have applied the CSL theory to inflation. Since 
the CSL scenario addresses the measurement problem in 
quantum mechanics, it is a priori relevant to explain how 
the wave-packet reduction took place in the early uni- 
verse, in the absence of any observer. Assuming that 
the wavefunction has to collapse on an eigenstate of 
the Mukhanov-Sasaki operator, we have computed the 
scalar power spectrum of cosmological perturbations and 
studied the dynamics of the wavefunction collapse. We 
have found that, in order to preserve the scale invari- 
ance of the power spectrum, it is necessary to fine-tune 
the parameter 7 which controls the amplitude of the 
CSL corrections. Typically, depending on which tempo- 
ral gauge is chosen (see the appendix A), we have found 
that the dimensionless parameter that can be constructed 
out of 7 must be smaller than exp (—a few AA^*), where 
ANt ~ 50 — 60 is the number of e-folds spent by the rel- 
evant modes outside the Hubble radius during inflation. 
We have also found that the time available during the 
inflationary phase is sufficient in order for the pertur- 
bations wavefunction to collapse. However, due to the 
smallness of 7, the spread of the final wavefunction is 
too important, rendering the collapse process not suffi- 
ciently efficient. Therefore, under the assumptions made 
in this paper, it seems fair to claim that the collapse theo- 
ries cannot solve the inflationary "macro-objectiflcation" 
question. 

The conclusions drawn above may not be as drastic 
as they appear at flrst sight, because they are subject 
to some assumptions, and in particular the choice of 
the "collapse operator" as the Fourier space Mukhanov- 
Sasaki Vk variable: all cosmological predictions made to 
date are based on this variable, rendering this choice very 
sensible, but it is by no means unique (see, e.g. the dis- 
cussion in Sec. VA). Moreover, can be understood 
as a quantum fleld living in a curved spacetime, so it 
should be treated by a quantum field theory version of 
the CSL mechanism. The present state-of-the-art of this 
subject technically forbids such a direct treatment, hence 
our simplifying hypothesis. Could it be that a full rela- 
tivistic version of CSL, reproducing the many successes of 
quantum field theory and of the ensuing particle physics, 
is needed before we can even embark in examining cos- 



mological perturbations? We doubt so, because cosmol- 
ogy, contrary to ordinary quantum field theory, is en- 
dowed with a preferred timelike direction that renders 
the "time-dependent Minkowski approximation" accu- 
rate enough for all practical purposes. It is left for future 
investigations to verify that the potential problems raised 
and stringent constraints obtained in this work could be 
naturally solved in a more general, yet-unknown, frame- 
work. 

There are other questions that could be the subject of 
further works. In particular, there is the issue that energy 
is not conserved in the CSL theory. In the case of the 
harmonic oscillator, this is expressed through Eq. (141). 
In the case of cosmological perturbations, it is easy to 
show that this leads to 
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The CSL contribution can easily be integrated and gives 
(^fc)lcsL — 7^ini/2 at the end of inflation. Expressed in 
terms of the Hamiltonian rather than the Hamiltonian 
density, one arrives at 



(H) 



fc^dfc, 



(242) 



which is infinite. It does not come as a surprise as it 
is known that the CSL Tomanaga-Schwinger equation 
precisely leads to this type of divergences [67, 08]. It 
could be regularized by introducing an ultraviolet cutoff 
although we notice on the above equation the weird prop- 
erty that the infinite integral is over comoving wavenum- 
bers rather than over physical ones. This energy non- 
conservation should cause a continuous increase of energy 
density during inflation. It is interesting to notice that 
it cannot occur at first order in the perturbations since 

fc) ) — 0- This means that it will be important at 
second order only. Then, it would be important to quan- 
tify this effect and, in particular, to compare it to the 
background energy density ~ H^M^^ in order to check 
whether this leads to a backreaction problem. 

Another point is that we have shown that the power 
spectrum, contrary to what happens in the standard case, 
remains a time-dependent quantity, i.e. still evolves with 
time on large scales during inflation. It is therefore not 
obvious that evaluated at the end of inflation is ex- 
actly the power spectrum that should be used at recom- 
bination. In fact, what happens just after the end of 
inflation is of great interest for the cosmological conse- 
quences of CSL. Indeed, just after inflation, the stages 
of pre-heating and re-heating begin [96-98]; this is also 
shown in Fig. 4. During this phase of evolution, the 
inflaton field oscillates at the bottom of its potential, 
(p{t) oc sin(mt -I- A)/{mt) where A is a phase and m the 
mass of the inflaton (in the case of power-law inflation, 
the potential has no minimum and, therefore, can only 
be used to describe the slow-roll regime. Here, we assume 
that the potential can be approximated by rn^ip^ in the 
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vicinity of the minimum). In this case, the equation of 
motion (12) for the Mukhanov-Sasaki variable takes the 
form of a Mathieu equation [99]. As is well known, this 
equation possesses unstable solutions when the parame- 
ters falls in the resonant bands. In the case of inflation, 
one can show that the large-scale perturbations are in 
the first instability band which makes growing and 
staying constant [99, 100]. In the CSL case, the corre- 
sponding Mathieu equation would read 



d^Ufe 



[Ak ~ 2(7Cos(2z + 2A)] Vk = 0, 



(243) 



where z = mt -\- 7r/4, Ce, te denoting the scale factor and 
the cosmic time at the end of inflation and with 
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Since g ^ 1, in the regular case when 7 = 0, the condi- 
tion to be in the first resonant bands, 1 — g < Ak < 1 -l-g, 
is equivalent to < k/a < ^SHm. In the CSL case, 
the coefficient Ak becomes complex. Therefore, in order 
to determine the corresponding Floquet index, it now 
becomes necessary to study the instability chart of the 
Mathieu equation in the complex domain. Although this 
is beyond the scope of this paper, this is certainly a sub- 
ject worth investigating. In particular, it would be inter- 
esting to see whether the instability is enhanced in this 
case as one can, maybe naively, suspect. If so, maybe the 
preheating stage can put even more stringent constraints 
on the parameter 7. 

We have seen that the study of the CSL cosmological 
perturbations is in fact equivalent to the study of the CSL 
parametric oscillator (i.e. an harmonic oscillator with a 
time-dependent frequency) . The previous discussion sug- 
gests that it would be interesting to investigate the case 
of a parametric oscillator in the presence of a resonance 
in the CSL framework. In quantum field theory, this is a 
common situation and typical examples are the dynam- 
ical Schwinger effect [115] (the analogy between cosmo- 
logical perturbations and the Schwinger effect was dis- 
cussed in Ref. [15]) or the dynamical Casimir effect [116] 
which was recently observed for the first time [117] in 
the laboratory. In fact, if we want to avoid the objection 
that the quantum field CSL theory is not yet ready, it 
would be even more interesting to find a non-relativistic 
system governed by a Mathieu equation and to investi- 
gate its behavior within the CSL theory. We believe that 
all the equations presented in the present article can be 
straightforwardly applied to this case. Here, we suggest 
that a Paul trap [IIX] could be such an example. As 
for the inflationary preheating, we expect the coefficients 
of the Mathieu equation to become complex because of 
the — 2«7 term. This will probably make the system ex- 
tremely unstable and, as a consequence, it will probably 
be possible to put very tight constraints on the value of 
7. We hope that this case will be treated in details soon. 
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Appendix A: "Gauge Invariance" of the CSL Power 
Spectrum 

In section V A, we discussed the choice of the "collapse 
operator" , i.e. the operator that appears in the non-linear 
and stochastic part of the modified Schrodinger equa- 
tion. In principle, this operator should be determined by 
a more fundamental theory. However, the CSL model is 
just a phenomenological approach and the "collapse op- 
erator" is just put by hand in order to match what we ob- 
serve when an experiment or an observation is performed 
(the position of a spot in a detector, the energy density of 
a field, etc.). In the case of the cosmological primordial 
perturbations, we have argued that the Mukhanov-Sasaki 
variables Vk is the most sensible choice. But this variable 
often appears factorized by a background quantity, typi- 
cally a power of the scale factor 0(77). Therefore, instead 
of Vk, one could very well choose the collapse operator to 
be h (a) Vk, where his a priori an arbitrary function of the 
scale factor a. After all, Vk and h (a) Vk share the same 
eigenspectrum and drive the system towards the same 
target states with the same probabilities. But the point 
is that, a priori and as is discussed in detail below, this 
does not lead to the same solution for the mode function 
fkiv) and, therefore, a priori, for the power spectrum. 

In fact, this issue is related to an even more funda- 
mental problem. Indeed, one could claim that the con- 
formal time r] used in this paper to write the modified 
Schrodinger equation is not the physical one and that 
one should use instead, say, the cosmic time t (of course, 
the discussion also applies to any other time variables re- 
lated to T] through a transformation that depends only on 
the background). In fact, a choice of time is equivalent 
to a choice of h since it has the same effect on the mod- 
ified Schrodinger equation. And, of course, as already 
mentioned, one could worry that different choices lead 
to different predictions. Therefore, the phenomenologi- 
cal approach used in this article suffers from what can 
be called a "temporal gauge" problem. This problem 
probably originates from the fact that the CSL equation 
is not covariant under diffeomorphisms (contrary to the 
standard theory of cosmological perturbations). 

In this appendix, we investigate this question, showing 
the remarkable property that the conclusions obtained 
in this paper for h (a) = 1 are in fact valid for any other 
functions h. It is true that the detailed shape of the power 
spectrum depends on the gauge but its global properties 
are independent of the choice of h. This means that, a 
priori for any h allowing meaningful initial conditions, the 
power spectrum of cosmological perturbations has a bro- 
ken power-law shape, with = 1 at small wavelengths 
and rig = 4 at large wavelengths. As a consequence, the 
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requirement of moving the non-scale-invariant part of the 
spectrum beyond the Hubble radius today always leads 
to extreme constraints on the parameter 7. 

Let us now consider the modified Schrodinger equation 
of motion for in the CSL picture, with spontaneous 
localization on the h (a) Vk eigenmanifolds. It reads 



d*« 



(Al) 



and a similar equation for This equation should be 
compared with Eq. (142), the only difference being that 
the operator is now multiplied by h{a). Parameteriz- 
ing vj/fc as in Eq. (143) using again fife = ^if'k/ (2./fe), one 
is led to the following equation for the mode function 

[^2(77,fe)-2z7/i2(a)]/fe = 0. (A2) 

This expression should be compared with Eq. (152): as 
expected, the only difference is that an extra /i^ (a) ap- 
pears in front of the 7 term. For simplicity, let us choose 
/i to be a simple power law and let us assume the in- 
flationary dynamics to be close to a de Sitter Universe 
a{ri) ~ —lo/rj. Then, the mode function can be re- 
expressed as 



77^ 



2«7a^ ./fc = . 



(A3) 



If p < 0, the Bunch-Davies vacuum state cannot be 
chosen at the onset of inflation since the fc^ term does 
not dominate in the parenthesis. This means that one 
must work with p > 0. In this paper the case p ~ Q 



[i.e. h (a) = 1] has been studied, hence one only needs 
to study the p > cases. It is interesting first to notice 
that the cases p> provide a natural amplification phe- 
nomenon depending on the physical length of the mode 
since the amplitude of the term proportional to 7 now 
increases as the mode is stretched by the growth of the 
scale factor. This is consistent with the physical intuition 
which tells us that the collapse should occur for macro 
extended objects only, li p > 2, the term proportional to 
7 dominates the dynamics at the end of inflation, when 
kr] goes to 0, and one can expect the power spectrum 
scale invariance to be destroyed. Therefore, if p is an in- 
teger, we are left with the cases p = 1 and p — 2 that we 
now study. 

li p — 1, the general solutions of Eq. (A3) can be ex- 
pressed in terms of Whittaker functions W^^^ (z) [101, 
102] as 

fkiv) = C!kW^eo/k,3/2 {2ikT]) + DkW_^eg/k,3/2 {-^ikif} , 

(A4) 

where Ck and Dk are integration constants that can 
be determined by choosing the Bunch-Davies vacuum 
state for the initial conditions. This leads to Ck = 0. 
Then, in the limit where fcry goes to 0, '^eVLk{vi) can 
be Taylor expanded, and this provides a simple ex- 
pression for this quantity. In particular, we find that 
5Rer2fe/fc — 7^o/(2fc) + C'(^??), showing that, in this case, 
the spread does not divergence in the large-scale limit and 
that, as a consequence, the localization of the wavefunc- 
tion becomes much more accurate. Moreover, since the 
inverse of SRe fife is basically Vc^ , this allows us to calculate 
the power spectrum, provided we push the expansion to 
higher orders. One obtains 



V^ik) 



47 
k 



kf) \ k 
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loe 



fco 



A TV. 



^C(^)|stand' 

m. 



where stand standard power spectrum (47), and 

where g {x) is defined by 

— = l+2,x~?,x^~x^-2x (1 - a;2) U {2 + x) - 2iP (1)1 , 
g{x) 

(A6) 

i/j (x) being the digamma Euler function [101, 102]. Let 
us notice that, in Eq. (A5), we have sometimes intro- 
duced the quantity £o"f/k. Of course, the most conve- 
nient way of dealing with this quantity is to express it as 
{£Qj/ko)ko/k such that the dimensionless small param- 
eter ioj/ko explicitly appears. The spectrum given by 
Eq. (A5) should be compared with the one obtained in 
Eq. (171) with the choice h = 1. They share the same 
broken power-law structure, with a scale-invariant part 



rig ~ 1 at small scales and a branch with Ug = 4 on large 
scales. This spectrum is displayed in Fig. 6 for different 
values of the parameter £07/^0- 

The break in the power spectrum appears at fc^/fcp ~ 
^07/^06^^^*. Therefore, in order for the non-scale in- 
variant part of the power spectrum to be outside the 
Hubble radius, one must have 

74 ^ ^-2AN, ^ ^0-53 (AT) 

ko 

This equation should be compared to Eq. (177). We see 
that, in the present case, we also obtain a constraint 
that can be considered as "extreme". In other words, it 
seems that a very important fine-tuning is necessary to 
maintain the consistency of the CSL predictions with the 
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CMB observations. We also notice that, instead of 7/A:q, 
it is now the combination "fio/ko that is constrained. 
Of course, this is just the consequence of the fact that, 
as aheady discussed, changing the coUapse operator can 
change the dimension of the parameter 7. In some sense, 
we face again the discussion of the temporal gauge issue. 

Let us now turn to the case p — 2 in Eq. (A3). The 
general solutions of this equation can be expressed in 
terms of Bessel functions with a complex order [101, 102], 
namely 

i-kv) , (A8) 



where and are integration constants that can 
be determined by requiring, as usual, the initial state 
to be the Bunch-Davies vacuum. This leads to Cfc = 
-Dfce^^'^/Vi+s/s^T^o. In the hmit where kr] goes to 0, 
SRe fife can be Taylor expanded and, at first order in the 
parameter j£q, the power spectrum reads 



Vcik) 



27T 
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V^k) 



stand ' 



(A9) 



The formula (A9) should be compared with Eqs. (171) shape, with a scale invariant part on small scales and 
and (A5). Again, the power spectrum has the same a non-invariant branch with rig = 4 on large scales. 
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This is clearly seen in Fig. 7, where the spectrum (A9) 
is represented for different values of the parameter j£q. 
The break in the power spectrum appears at k^/k^ ~ 
7^g/3e'^^^*. Therefore, in order for the non-scale in- 
variant part of the power spectrum to be outside the 
observational window, one must require that 

j£l « e-''^'^* ^ 10-^^ . (AlO) 

Again, we can consider the above constraint as a fine- 
tuning. It is also interesting to notice that, contrary to 
Eqs. (177) or (A7), Eq. (AlO) involves physical quantities 
only. This is because, when p = 2, the CSL correction 
that should be compared to the comoving wavenumber 
squared is oc 70^, see Eq. (A3). In other words, 7 should 
now be compared to the physical wavenumber. If we take 
£q ~ lO'^^pj, which comes from the CMB normalization, 
then one arrives at 7 ^ 10^*^. 

Let us conclude this appendix by noticing that the 
above results are in fact generic and do not depend on 
the value of p. Technically, the power spectrum is ob- 
tained by taking the super-Hubble limit of the mode func- 



tion fk (r?) , by inserting it in the expression of Sff e fife = 
SRe [— «/fc/(2/fc)] and by retaining only the leading order 
in krj. In the standard case, the leading terms of the mode 
function expansion turn out to cancel out in 5fte Qk, leav- 
ing an expression which precisely gives a scale-invariant 
power spectrum. This cancellation originates from the 
fact that the Wronskian is conserved. In the CSL case, 
the fact that 7 7^ implies that this symmetry no longer 
exists, and, as a consequence, the nice cancellations men- 
tioned above no longer show up and scale invariance is 
immediately broken. In some sense, the fact that the 7 
term destroys the scale invariance of the power spectrum 
does not come from the fact that its presence modifies 
the time dependence of the effective frequency (the value 
of p or the choice of h), but is rather due to the fact that 
it makes the effective frequency a complex quantity. We 
conclude that modifying the definition of the "collapse 
operator" by multiplying it with a background function, 
despite changing the dimension of 7, always constrains 
this parameter to be extremely fine-tuned. 
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